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FOREWORD 


This  is  the  Final  Report  under  Contract  F19628-79-C-0033  with  the 
Meteorology  Division,  Air  Force  Geophysics  Laboratory.  The  contract 
extended  from  December  1,  1978  to  November  30,  1980. 

The  report  consists  of  narrative  descriptions  of  work  accomplished 
during  those  two  years  and  two  technical  papers  reporting  research 
achievements.  Work  which  was  described  or  reported  on  previously 
(in  Scientific  Reports  Nos.  1,  2,  and  3)  is  referenced  in  the  appro¬ 
priate  places  but  is  not  repeated  except  in  brief  summary  form. 

The  contract  called  for  effort  in  five  discrete  areas,  reported 
upon  in  the  monthly  R&D  Status  Reports  for  the  last  two  years  under 
these  headings: 

Numerical  Weather  Prediction 
Automated  Weather  System 
Satellite  Meteorology 
Me I DAS  Operation  and  Enhancement 
Weather  Radar  Techniques. 

This  Final  Report  is  organized  similarly. 

Personnel  associated  with  these  areas  at  any  time  during  the  contract 

are: 

Numerical  Weather  Prediction  -  Donald  C.  Norquist,  M.S. 

Peter  J.  Mansfield,  Ph.D. 

Isidore  M.  Halberstam,  Ph.D. 

Gail  M.  Bertolini,  A . B . 

Dorothy  C.  Moran,  A.B. 

Scott  D.  Hamilton 


Automated  Weather  System 


Satellite  Meteorology 


McIDAS  Operation  and 
Enhancement 


Joan  M.  Ward,  A.B. 

Russell  C.  Dengel,  B.S. 
Maureen  F.  Hampton,  B.A. 

Gail  A.  Dengel,  M.S. 

Robert  P.  d'Entremont,  B.S. 
Michael  E.  Niedzielski,  B.S. 
Lucy  K.  Aharonian,  A.M. 
Laurie  E.  Engels 
John  M.  Powers,  B.A. 

Barry  A.  Mareiro 
Judith  E.  Larkin 
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Weather  Radar  Techniques  -  Alan  R.  Bohne,  Ph.D. 

John  D.  Tuttle,  M.S. 

Wayne  S.  Collamore,  B.S. 

Chen-Hung  Chang,  M.S. 

Glenn  R.  Smythe,  M.S. 

Principal  Investigator  was  Alan  M.  Gerlach,  Ph.D. 

Subcontracts  with  Harris  Corporation  Computer  Systems  Division  and 
Perkin-Elmer  Customer  Service  Division  provided  on-site  maintenance  ser¬ 
vices  for  the  McIDAS  and  Weather  Radar  Data  Processing  System  respectively. 

Reports  were  prepared  by  the  scientists,  engineers,  and  mathemati¬ 
cians  identified  in  the  Table  of  Contents  with  their  technical  reports. 
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A.  THE  EFFECT  ON  FORECAST  ACCURACY  OF  VARYING  SPATIAL  RESOLUTION  IN  A 
LIMITED-AREA  NUMERICAL  WEATHER  PREDICTION  MODEL 

I.  Purpose  of  Study  and  Model  Equations 

One  of  the  many  challenges  facing  scientists  engaged  in  numerical 
weather  prediction  modeling  is  that  of  minimizing  prediction  error  while 
maximizing  computational  speed  and  efficiency.  In  finite  difference  models 
such  as  the  six-layer  Limited-Area,  Fine-Mesh  (LFM)  model  initially  formu¬ 
lated  by  Shuman  and  Hovermale  (1968)'*",  the  large  data  lat  tice  on  which  the 
atmospheric  primitive  equations  are  solved  numerically  requires  a  large 
computer  storage  capacity.  In  addition,  multitudinous  calculations  are 
made  at  each  time  step  of  the  model  simulation  so  that  relatively  large 
processing  times  are  involved.  Although  a  finer  grid  resolution  would 
yield  better  approximations  to  the  derivatives  in  the  primitive  equations, 
it  would  also  necessitate  a  proportionately  greater  outlay  of  computational 
storage  and  processing  time.  Furthermore,  the  necessary  spatial  interpola¬ 
tion  of  the  initial  data  and  boundary  values  onto  the  finer  grid  could 
significantly  alter  the  input  data  and  degrade  the  resulting  forecast. 

This  Section  describes  an  effort  to  investigate  the  effect  of  varying 
the  spatial  resolution  on  forecast  accuracy  by  using  a  proven  operational 

forecast  model.  The  National  Meteorological  Center  (NMO)  six-layer  LFM 

2 

model  as  documented  by  Gerrity  (1977)  was  coded  for  the  Air  Force  Geo¬ 
physics  Laboratory  (AFGL)  computer.  The  model  was  formulated  for  computa¬ 
tions  on  a  33  x  33  point  horizontal  grid  at  each  of  six  prognostic  levels 
in  the  vertical,  then  was  modified  to  perform  the  same  computations  on 
grids  containing  twice  and  four  times  the  number  of  grid  points  at  each 
level.  Analyzed  data  at  approximately  the  horizontal  density  of  the  33 
x  3)  grid  were  interpolated  to  provide  finer  resolution  initial  and  bound- 

^Shuman,  F.G.,  and  J.B.  Hovermale,  1968:  An  operational  six-layer  primitive 
equation  model,  J.  Applied  Meteor.,  7,  525-547. 

^Gerrity,  J.F. ,  Jr.,  1977:  The  LFM  Model  -  1976:  A  documentation .  National 
Meteorological  Center,  NOAA  Technical  Memorandum  NWS,  NMC  60. 
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ary  values  for  the  modified  model.  The  goal  of  this  study  was  to  deter¬ 
mine  the  optimal  spatial  resolution  for  a  given  input  data  density,  in 
order  to  achieve  maximum  accuracy  at  minimum  computational  time  and  ex¬ 
pense. 

A  f  ill  discussion  of  the  method  used  in  the  model  to  solve  the 
primitive  equations  was  included  in  an  earlier  report  by  the  author 
(Gerlach,  ed. ,  1980)  \  The  primitive  equations  in  Cartes. an  form  for  a 
conformal  projection  of  the  earth  are  repeated  below;  the  symbols  are  ex¬ 
plained  in  the  earlier  report: 
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(1) 


(2) 

(3) 

(4) 


(5) 

(6) 


(7) 

(8) 


The  coarsest  resolution  grid  (33  x  33  points  on  each  of  six  prog¬ 
nostic  levels  and  one  computational  level)  was  set  up  to  cover  that  seg- 


3 

Gerlach,  A.M.,  ed.,  1980:  Computer-based  weather  research.  AFGL-TR-80- 
00G9.  Scientific  Report  No.  2_,  Contract  F19628-79-C-0033 ,  pp.  5  -  19. 
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ment  of  the  NMC  octagonal  qrid  which  completely  includes  the  contiguous 
United  States  (see  Figure  A-l).  The  finer  resolutions  of  the  model  in¬ 
volved  grid  points  within  this  same  region  but  at  twice  and  four  times  the 
density  at  equally-spaced  intervals.  Finite  differences  were  calculated 
using  the  method  described  by  Shuman  and  Hovermale  (1966) .  The  input  data 
were  prepared  so  as  to  accommodate  both  second  and  fourth  order  spatial 
finite  differences,  but  only  the  second  order  scheme  was  actually  coded  in 
the  model.  The  model  used  a  center-difference  in  time  (leapfrog)  method 
except  at  the  first  time  step  where  forward  difference  was  used. 

II.  Implementation  of  Variable  Resolution  in  the  Model 

The  computer  code  for  the  33  x  33  x  7  grid  point  model  described  in 
this  report  was  formulated  following  a  procedure  identical  to  that  des¬ 
cribed  in  the  Gerrity  documentation.  The  only  exception  is  the  absence 
of  the  time-smoothing  technique  used  in  the  NMC  model.  This  was  not  in¬ 
cluded  because  it  was  thought  that  this  effect  might  mask  the  effect  of 
changing  the  model  resolution.  After  successful  testinc  in  this  form,  the 
model  code  was  modified  to  permit  the  model  computations  to  be  carried  out 
on  65  x  65  x  7  and  129  x  129  x  7  grid  lattices.  Since  the  originally 
formulated  coarse  resolution  version  of  the  model  required  essentially  all 
of  the  core  memory  available  on  the  AFGL  CDC  6600  computer,  a  simple  in¬ 
crease  of  storage  array  sizes  to  accommodate  the  higher  resolutions  was 
not  possible.  Instead,  the  code  was  redesigned  to  carry  out  the  computa¬ 
tions  for  the  larger  grids  in  35  x  35  x  7  grid  point  sectors.  Four  such 
sectors  were  required  for  the  medium  resolution,  and  16  were  needed  for 
the  fine  resolution.  In  this  way,  the  core  memory  requirement  remains 
constant  with  change  of  model  resolution,  and  only  the  computation  time 
increases  accordingly. 

In  the  code.  Equations  (1)  -  (5)  are  solved  numerically  in  sequence 
for  each  sector,  after  which  the  sector's  results  are  sent  to  mass  stor¬ 
age  before  the  computations  for  the  next  sector  begin.  The  sectors  are 
processed  in  a  distinct  sequential  order,  beginning  with  the  southwestern- 
most,  proceeding  northward  to  the  northwesternmost,  then  to  the  sector 
just  east  of  the  initial  sector  and  northward  again,  and  coi  j.nuing  with 
this  process  until  the  northeasternmost  sector  is  completed. 
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Fiquro  A- I .  The;  NMC  octaqona J  qrid.  Square  outlined  reyion  is  the 
limited  ai  ea  considered  in  this  study.  Coarsest  resolution  considered 
is  I  ..ice  that  of  thi:  qrid.  (Adapted  from  Shuman  and  Hovermale,  1968.) 
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After  all  newly  predicted  internal  values  for  the  stctors  are  stored 
in  mass  storaqe,  external  boundary  values  for  that  time  step  are  included 
with  the  computed  values  and  the  entire  sectors'  values  ire  stored.  Then, 
sector  by  sector,  the  algorithms  for  precipitation  proce;  ses  are  imple¬ 
mented  in  which  all  values  of  velocity,  potential  temper,  ture,  and  pre- 
cipitable  water  in  the  sectors  are  subject  to  possible  adjustment. 

Finally,  the  check  for  super  adiabatic  lapse  rate  is  carried  ou:  and  the 
"dry  adiabatic  convection  adjustment"  is  implemented  for  each  column  in 
which  such  instability  is  found.  It  was  necessary  to  limit  the  number  of 
attempts  to  remove  instability  completely  from  a  particular  column}  the 
limit  was  placed  at  300  tries,  in  this  study.  Upon  completion  of  the  super 

adiabatic  check,  the  newly  predicted  and  adjusted  values  of  — ,  -,  0, 

‘  '  m  m  JO 

and  W  are  again  sent  to  storage  and  the  precipitation  adjustment  proce¬ 
dure  is  begun  for  the  next  sector. 

In  order  to  perform  the  finite-difference  operations  on  the  rows  or 
columns  which  form  the  border  between  two  adjoining  sectors,  three  rows  or 
columns  common  to  both  sectors  are  involved  in  the  computations  for  both 
sectors.  This  overlap  is  necessary  to  provide  a  boundary  row  in  order  to 
compute  internal  values  for  each  sector.  What  is  used  as  a  boundary  row 
for  computations  in  one  sector  is  the  last  row  for  which  computations  were 
made  in  the  adjoining  sector.  The  extra  storage  space  required  for  these 
rows  was  provided  by  having  originally  coded  the  model  for  a  35  x  35  grid 
which  when  doubled  or  quadrupled  allows  space  for  up  to  70  x  70  or  140  x 
140  respectively,  enough  to  accommodate  even  a  fourth  order  differencing 
scheme  if  it  were  to  be  incorporated  into  the  model  in  the  future. 

III.  Preparation  of  Initial  and  Boundary  Values 

The  model  formulation  followed  exactly  the  methodology  used  in  the 
NMC  model.  This  was  not  true  for  the  preparation  of  the  input  data  for 
the  model,  however.  The  differences  are  due  largely  to  the  nature  and 
spatial  distribution  of  the  atmospheric  data  chosen  as  input  for  the 
model.  The  data  for  initial  values  and  boundary  values  are  taken  from 

4 

the  FGGE  Level  III-A  data  set  .  The  observed  data  had  been  objectively 
analyzed  to  produce  these  data  values  located  at  2.5°  x  2.5°  latitude- 

4 

Obtained  from  World  Data  Center-A,  Na'  onal  Climatic  Center,  Ashe¬ 
ville,  NC  28801. 
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longitude  intervals  and  at  standard  mandatory  pressure  levels.  Values 
for  zonal  and  meridional  wind  components,  temperature,  geopotential  height 
up  to  50  mb,  and  relative  humidity  up  to  300  mb  are  available  at  each  grid 
point.  In  addition,  surface  values  of  temperature  and  geopotential  height 
are  included  in  the  data  set,  along  with  mean  sea  level  pressure  and  tropo- 
pause  temperature  and  pressure.  Surface  albedo,  drag  coefficient,  terrain 
elevation,  and  sea  surface  temperature  values  are  also  provided  and  are 
used  Ln  the  parameterization;:  included  in  the  model.  Since  data  values 
from  the  end  time  of  the  forecast  are  used  to  compute  boundary  values, 
both  initial  and  final  data  values  are  processed  through  all  steps  leading 
up  to  preparation  of  the  boundary  values. 

(A) .  Vertical  Interpolation  from  Pressure  Levels  to  0  Levels 

Once  the  desired  data  are  extracted  from  the  larger  data  set,  the 
first  step  in  preparing  the  data  for  model  input  is  the  vertical  inter¬ 
polation  of  data  values  from  mandatory  pressure  levels  to  a  levels  re¬ 
quired  by  the  model.  The  code  which  performs  this  interpolation  pro¬ 
cesses  each  vertical  column  separately.  The  latitude  and  longitude  of 
the  column  are  used  along  with  the  velocity  components  to  obtain  the  velo¬ 
city  com|onents  along  grid  lines. 

At  this  point  in  the  code  the  pressure  at  each  of  the  0  levels  is 
calculated  in  order  to  obtain  for  each  of  the  three  domains  (tropo¬ 
sphere,  stratosphere,  and  computational  layer)  according  to  the  definition 
of  ij  in  each  domain  (see  Gerlach,  ed. ,  1980) .  Although  mean  sea  level 
pressure  is  provided  for  each  column,  it  is  not  used  to  obtain  surface 
pressure  where  the  terrain  lies  above  sea  level.  Attempts  to  extrapolate 

a  value  for  p  _  in  such  cases  using  p  ,  in  the  hydrostatic  equation  of- 
sfc  msl 

ten  yielded  a  surface  potential  temperature  greater  than  that  for  the  a 
level  immediately  above  it.  As  a  result,  surface  pressure  is  determined 
usinq  the  hydrostatic  assumption  in  the  form 

p  =  p  exp  f( $  -  3>  ,  )/Rt]  (9) 

whore  T  is  the  arithmetic  average  of  the  surface  temperature  and  the  tem¬ 
perature  at  the  first  mandatory  pressure  level  (p  )  above  the  surface. 

m 

The  surface  pressure  (p^)  and  the  given  tropopause  pressure  (p^)  are  then 
used  to  generate  pressure  values  p^ ,  p^ ,  and  p^  for  the  intermediate  0 
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levels  as  defined  in  the  previous  report  (Gerlach,  ed.  ,  1980).  The 

corresponding  values  of  the  Exner  function  at  each  of  these  levels  is 

then  obtained  using  Equation  (7)  . 

To  calculate  the  value  of  geopotential  at  the  tropopause  (  4  )  the 

5 

geopotential  is  assumed  to  be  a  quadratic  function  of  tt  in  the  strato¬ 
sphere: 


4(TT)  =  <J>  -  A(TI  -  7T  )  -  71 

5  5  2 


V 


(10) 


Differentiating  (10)  with  respect  to  tt,  it  can  be  seen  that  A  =  0^  and 

B  =  c  [0  -  0(TT)1  /  (tt  -  tt)  .  Eq.  (10)  is  solved  for  4*  using  data  from  the 

poo  5 

next  mandatory  level  above  the  tropopause,  or  in  the  case  of  <  lOOmb, 
data  at  100  mb  are  used  to  obtain  4^.. 

The  value  of  geopotential  at  the  top  of  the  stratospheric  domain 
(4^)  is  obtained  using  the  same  assumption.  Here  the  quadratic  assump¬ 
tion  takes  the  form 


4?  =  4  -  A(T7?  -  TT)  -  |(7T7  -  TT)2  (11) 

where  TT,  4  refer  to  the  values  of  the  respective  variables  at  the  mid-level 

(in  it)  of  a  layer  bounded  by  mandatory  levels  at  or  above  the  tropopause. 

Using  "U"  and  "L"  to  denote  the  known  values  at  the  top  and  bottom  of 

this  layer  respectively,  we  have  TT  =  ^(tt  +  7T  )  ,  A  =  -  (4  -4  )/fn  -  tt  )  , 

Li  U  LULU 

and  B  =  c  (0.  -  0,,) / (tt  —  TT  )  .  The  expression  for  4  can  be  derived  as 
p  L  U  __L  _  u 

follows:  substitute  4,  tt  for  4,tt  in  Equation  (10)  ,  and  replace  the  "5" 

subscript  in  the  result  with  "L".  Then  noting  from  the  definition  of  n 

that  tt  —  tt  =  *5(tt  —  tt  )  ,  substitute  this  into  the  result.  Form  another 

L  U  L» 

equation  by  substituting  4  ,  tt  for  4,tt  and  4  ,  tt  for  4  ,  tt  in  Equation 

U  U  Li  L<  D  D 

(10) ,  and  combine  this  result  with  the  equation  just  obtained  bv  eliminat¬ 
ing  the  first  order  term  between  the  two  equations.  The  resulting  ex¬ 
pression  for  4  is 

4  =  *s(4t  +  4  )  +  |(ttm  -  tt)2  .  (12) 

In  the  current  study  the  "U"  subscript  denotes  values  of  the  subscripted 
variables  at  100  mb  and  "L"  denotes  values  at  150  mb  when  the  tropopause 
pressure  is  greater  than  150  mb.  Values  at  100  mb  are  used  for  the  vari¬ 
ables  with  "U"  subscripts  and  tropopause  values  for  "L"  subscripted  vari¬ 
ables  when  p_  <  150  mb. 

D 

The  symbol  4TOp  represents  the  geopotential  at  the  top  of  the  model 
atmosphere,  a  value  which  is  chosen  arbitrarily.  Thus,  4^  =  4 
-c  6^7*  w^ere  t^ie  constant  potential  temperature  in  the  oomputa- 


tional  layer.  To  derive  an  expression  lor  n?  in  terms  of  known  quanti¬ 
ties,  this  is  rewritten  in  the  form 

*7  '  *TOP  -  CP°c’  -  -  «  '  (U> 

Subtracting  (11)  from  (11)  yields 

(tt7  "  ¥)2  -  [2(cp0c  -  A)/b](t,7  -  •'„')  -i  2C/B  =  0  (14) 

where  C  =  <J>TOp  “  4>  "  c  6cii  -  Solving  for  using  the  quadratic  formula 

and  choosing  the  negative  sign  gives 

1T7  =  IT  +  - -  =  •  (15) 

ce  -  A  +  v ( c  G  -  A)  -  2BC 

P  C  p  c 

The  first  step  taken  in  solving  (15)  for  it  is  to  calculate  an 

arbitrary  value  for  0^  that  would  be  satisfactory  for  the  soundings  of 

all  grid  points  at  the  time  of  interest.  The  desired  value  of  0  is 

c 

chosen  using  the  sounding  data  from  the  southwesternmost  grid  point  in 

order  to  obtain  the  largest  value  of  s  .4249630111  (50  mb) .  The 

quantity  0^  =  (^TOp  ~  $) /c  it  has  its  lowest  values  in  the  tropics,  and 

the  lower  values  will  help  insure  real-valued  solutions  in  (15)  since  B 

is  always  positive.  Then  by  holding  0  constant  at  this  value  for  all 

c 

subsequent  soundings,  tt7  will  be  prevented  from  becoming  too  small  for 

other  tropical  soundings  or  too  large  (tt7  is  restricted  to  be  less  than 

.5180152937,  or  100  mb)  for  high  latitude  soundings. 

The  following  method  is  used  to  solve  (15):  starting  with  the  south- 

5  2  2 

westernmost  sounding,  choose  a  value  for  $TOp  equal  to  4.41  x  10  m  /sec 

(45  km),  calculate  0c  for  that  sounding,  round  the  result  up  to  the  next 

higher  value  that  is  evenly  divisible  by  ten  and  calculate  tt  .  If  the 

2  2 

discriminant  or  tt 7  become  negative,  increase  $TOp  by  9,800  m  /sec  (1  km) 

and  recalculate  tt_.  If  tt_  £  .4249630111  (50  mb),  save  this  value  of  0 
7  7  c 

for  all  soundings.  If  not,  increase  0  by  10  K  and  recalculate  tt_.  For 

c  7 

subsequent  soundings,  calculate  tt7  using  the  value  of  0^  obtained  for  the 

first  sounding.  If  the  discriminant  or  u_,  become  negative,  increase  $ _ 

2  2  7  ’  TOP 

by  980  m  /sec  (0.1  km)  and  try  again.  Once  a  real  positive  value  of 

is  calculated,  it  is  checked  to  see  that  tt  <  .5180152937  (100  mb).  If 

2  2  ' 

not,  lower  by  980  m  /sec  and  try  again.  Eventually  a  satisfactory 

value  of  tt7  will  be  obtained,  and  this  value  is  used  to  obtain  <J>7,  using 

Equation  (11),  as  well  as  to  obtain  tt  =  H(ttc  +  tt_)  . 

6  5  7 

Geopotential  at  all  other  diagnostic  levels  (k  =  2,3,4,  and  6)  is 
calculated  assuming  a  quadratic  relationship  in  the  natural  logarithm  of 
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pressure,  as  described  by  Gerrity  (1977;  see  Equations  A-2  through  A-9) . 

One  modification  was  made  to  the  scheme  for  calculating  $  :  when  p 

k=2  sfc 

is  less  than  the  next  mandatory  pressure  level  lying  below  Pj<_2'  t^ie 
surface  data  are  used  in  place  of  the  mandatory  level  data.  This  change 
combined  with  the  use  of  temperature  averaging  in  the  equation  for  p 

sfc 

eliminates  the  problem  of  9k=2  <  (calculated  at  a  later  step)  in 

mountainous  regions.  It  was  observed  that  the  calculation  for  ,  occa- 

k=6 

sionally  produced  a  value  for  6^  that  was  less  than  0r.  Therefore,  a 

D  5 

comparison  of  9,.  with  0_  was  included  in  the  code.  If  0  <  0r,  $  and 
6  5  6  5  TOP 

2  2 

<t>  are  both  increased  by  980  m  /sec  (to  increase  0  while  keeping  0_ 

7  6  7 

constant)  and  another  check  is  made.  When  this  instability  is  resolved, 
the  potential  temperatures  are  calculated  for  the  entire  column  using 
the  hydrostatic  equation,  except  6^  when  the  surface  pressure  is  less 
than  850  mb.  In  this  case,  9  is  obtained  using  a  linear  interpolation 
of  temperature  in  the  natural  logarithm  of  pressure  between  the  surface 
and  the  mandatory  level  above  it.  The  stability  of  the  column  is  checked, 
and  the  dry  adiabatic  convection  adjustment  described  by  Gerrity  is  used 
in  case  any  instabilities  occur.  Finally,  the  values  of  geopotential  are 
recalculated  from  the  stable  values  of  potential  temperature. 

The  final  step  in  the  vertical  interpolation  procedure  is  the 
generation  of  values  for  u,v  at  prognostic  sigma  levels  using  a  linear 
interpolation  in  7T.  This  is  done  by  putting  the  mandatory  pressure  level 
values  of  u,v  in  order  of  decreasing  tt,  finding  the  two  mandatory  levels 
which  bound  each  a  level  by  comparing  its  tt  value  with  the  mandatory 
level  values,  and  performing  a  linear  interpolation  between  the  two  to 
obtain  u  and  v  from  their  corresponding  values  at  the  mandatory  levels. 

The  u  and  v  values  at  70  mb  and  50  mb  are  used  to  interpolate  values  for 
all  a  levels  with  pressures  less  than  100  mb.  Values  for  the  specific 
humidity  q  at  a  levels  are  obtained  using  a  linear  interpolation  in  pres¬ 
sure  between  mandatory  levels. 

(B).  Interpolation  on  0  surfaces 

The  values  of  u,  v,  <J>,  and  q  for  the  appropriate  o  locations  at 

<30 

all  grid  points  are  passed  on  to  a  code  which  performs  a  series  of  inter¬ 
polations  at  each  0  location  for  each  variable.  In  determining  the  best 
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method  of  interpolating  values  onto  the  finer  resolution  grids,  bilinear, 
biquadratic,  and  trigonometric  interpolation  schemes  were  considered.  An 
earlier  investigation  of  three  methods  (trigonometric,  bicubic  spline, 
and  four-  point  averaging)  of  interjxslation  to  complete  a  staggered  grid 
(Norquis  ,  1979)  showed  the  trigonometric  method  to  be  somewhat  better 
than  the  other  two  in  preserving  the  spectral  components  of  the  original 
data  fie  d.  When  the  trigonometric  interfxalat ion  was  compared  with  the 
bilinear  and  biquadratic  methods  in  interpolating  onto  points  representing 
the  next  finer  resolution  grid,  no  similar  significant  advantage  was  ob¬ 
served.  Therefore,  since  the  bilinear  and  biquadratic  interpolation  al¬ 
gorithms  had  already  been  coded  for  the  finer  resolution  interpolations, 
they  were  used  in  the  current  study. 

Fir  ;t,  a  four-point  bilinear  interpolation  is  performed  in  which  the 
o  o 

data  values  at  each  2.5  x  2.5  latitude  grid  point  are  used  to  obtain 
values  on  a  35  x  35  grid  corresponding  to  the  grid  network  of  the  model. 
This  set  of  points  includes  both  the  33  x  33  grid  of  the  coarsest  reso¬ 
lution  plus  an  extra  row  of  grid  points  surrounding  it  which  are  used  for 
subsequent  interpolations  and  initialization.  After  this  interpolation, 
two  successive  biquadratic  interpolations  are  carried  out  to  obtain  the 
values  at  grid  points  representing  twice  and  four  times  the  resolution  of 
the  35  x  35  grid.  A  16-point  least-squares  biquadratic  interpolation  is 
used  in  which  a  slight  modification  is  required  when  a  grid  point  of 
higher  resolution  falls  on  a  grid  line  of  lower  resolution.  The  value 
of  such  u  point  is  obtained  by  averaging  the  two  least-squares  estimates 

on  the  two  sides  of  the  line.  After  all  interpolations  are  completed, 

Bp 

domain  values  of  and  the  values  of  <J>  at  diagnostic  levels  are  used  to 
evaluate  potential  temperature  0  at  the  six  prognostic  levels,  using 
Equation  (6) .  The  values  generated  for  the  computational  layer  are  aver¬ 
aged  and  this  average  value  is  assigned  to  all  grid  points  in  this  layer. 

Also,  precipitable  water  W  is  calculated  from  q  and  — •  and  reduced- 

da 

image  velocity  components  are  formed  by  dividing  the  velocity  component 
values  by  the  appropriate  grid  point  value  of  the  map  factor  m. 

5 

Norquist,  D.C. ,  1979:  An  investigation  of  three  methods  of  completing  a 
staggered  data  field  on  a  square  grid.  AFGL-TR-79-0197 .  Scientific  Re¬ 
port  No.  1,  Contract  No.  F19628-79-C-0033 . 
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(C) .  Wind  Initialization 


The  input  data  processing  continues  with  the  ini t iali zatio  i  of  the 

wind  field  at  each  prognostic  level  and  for  each  resolution.  1  lis  is 

done  by  assuming  that  the  given  reduced-image  wind  field  on  the  grid  can 

be  represented  by  a  vector  function  V  =  (— )i  +  (— ) "]  where  1  and  }  are 

m  m 

unit  vectors  parallel  to  grid  lines  and  positive  eastward  and  northward 
respectively.  It  is  further  assumed  that  V  can  be  considered  t  >e  sum  of 


V,  =  k  x  v<i> 

4> 


(16) 


-+d 


3x 


-*  i  ->■ 

where  k  is  the  unit  vector  in  the  vertical  pointing  upward,  V  =  i 
+  g'  the  horizontal  del  operator,  ijJ  is  the  scalar  stream  function,  and 

o  y 

X  the  scalar  velocity  potential  function.  Then  it  can  be  shown  that  the 
vorticity  £  is 

-4-  0 

(17) 


->  rt  _,2  , 

£=k‘VxV  =  Vl(t 
2 


where 


„2  -t  3 _  -f  3 _ 

V  =  1  — r  +  1  — r- 


.  ]  — —  is  the  horizontal  laplacian  operator,  while  the 
3x^  3y 


divergence  D  is  given  by 

D  E  ^  *  V  =  Vx  .  (18) 

In  order  to  formulate  a  wind  field  less  likely  to  contain  observation, 
analysis,  or  interpolation  errors  that  might  excite  spurious  motions  in 
the  forecast,  a  non-divergent  and  an  irrotational  wind  value  were  cal¬ 
culated  at  each  of  the  interior  grid  points  for  each  resolution  for  both 
the  initial  and  final  data  sets. 

The  non- divergent  component  values  were  obtained  by  first  calculating 

the  vorticity  L  =  -?—(—)  -  -*—(—)  numerically  for  each  grid  squaie  center 
dy  m  dx  m 

using  the  Shuman-Hovermale  (1968)  finite  difference  scheme.  The  next 
step  is  to  solve  (17)  numerically  to  obtain  values  of  ^  at  the  interior 
grid  square  centers.  Values  of  i|>  for  the  centers  of  the  exterior  row  of 
grid  squares  are  needed  to  do  this,  and  these  values  are  estimated  as 
follows . 

First,  assume  that  on  the  exterior  grid  square  centers  the  comjwnent 
of  the  wind  normal  to  the  grid  as  a  whole  is  due  only  to  the  nc n-d i vergent 


wind.  If  V  is  expressed  in  natural  coordinates  as 


*  --(I! 


(If  *  S)i 


where  s  and  n  are  unit  vectors  tangential  and  normal  to  the  boundary  of 

3y  -►  -> 

the  grid  as  a  whole,  then  this  assumption  means  that  ^  =  0  and  V  •  n 

34;  9n 

=  Lettinq  the  subscript  "i”  indicate  the  counter  for  the  exterior 

grid  squares,  we  have  in  finite  difference  form 


+  ^jf(As)  =  ijt  +  (V  •  n)  (As)  . 


The  normal  component  of  the  given  wind  at  the  sides  of  the  exter¬ 
ior  grid  squares  is  then  obtained  by  averaging  the  normal  component  values 
located  at  the  grid  points  on  both  the  interior  and  exterior  edges  of 
each  exterior  qrid  square  side.  Since  the  divergence  theorem  dictates 
that  the  sum  of  the  normal  components  of  the  wind  over  the  boundary  of 
the  region  equals  the  total  divergence  in  the  region,  this  sum  must  vanish 
for  non-divergent  flow.  This  sum  is  obtained  using  each  of  the  normal 
component  values  on  the  exterior  grid  square  sides  as  calculated  above, 
then  it  is  divided  by  the  number  of  such  normal  components  to  find  the 
value  each  normal  component  is  to  be  adjusted  by  (all  adjusted  equally) 
to  make  their  sum  vanish.  Then  Equation  (20)  is  used  with  these  adjusted 
values  to  obtain  estimates  of  in  the  exterior  grid  square  centers,  be¬ 
ginning  with  =  0  in  the  lower  left  corner  of  the  grid.  These 

boundary  values  of  ip  are  used  with  the  interior  values  of  £  to  solve  the 
Poisson  equation  (17)  numerically.  The  resulting  values  of  ip  in  the  in¬ 
terior  grid  square  centers  are  then  used  with  the  boundary  values  of  \p 

.  u  v  v  cM 

estimated  above  to  obtain  (— )  ,  =  -  and  (—)  .  =  -r1-  at  the  interior  grid 

m  if)  3y  m  p  3x 

points  using  finite  differences. 

The  calculation  of  the  irrotational  component  of  the  wind  at  interior 

3  u 

grid  points  is  done  in  an  analogous  way.  The  divergence  D  =  -r—  (— ) 

rv  dx  m 

3  v 

+  -^y(— )  is  calculated  for  each  grid  square  center.  In  obtaining  values 
of  x  for  the  exterior  grid  square  centers,  we  note  that  the  motion  is  ir- 
rotational  in  the  interior  only  if  the  circulation  around  the  region, 

9v 

ds  5  0,  this  condition  is  met 

ds 

if  the  line  integral  of  the  tangential  component  of  the  non-divergent  wind 
vanishes . 


given  by 


ds,  vanishes.  Since 
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In  the  procedure  used  in  this  study,  this  component  is  removed 
everywhere  on  the  boundary  in  the  following  way.  The  tangential  coi  po- 
nent  of  the  non-divergent  wind  on  the  grid  points  just  inside  the  e; terior 
grid  points  calculated  above  is  averaged  with  the  tangential  compon<  nt  of 
the  given  wind  on  each  exterior  grid  point  (in  which  it  is  assumed  that 
the  given  wind  is  largely  non-divergent)  to  obtain  an  estimate  of  V  •  s 
at  the  side  of  each  exterior  grid  square.  This  is  subtracted  from  the 
full  wind  value  at  this  same  location  obtained  as  in  the  non-di vergent 
wind  discussion,  to  obtain  V  •  s  at  each  point.  Then 


V  =  y  +  & 

xi+l  xi  As 


(As)  =  x-  +  (V  *  s) (As) 
■*  A 


<::i> 


is  used  to  obtain  values  of  x  at  each  exterior  grid  square  centers  begin¬ 
ning  with  =  ® '  T^e  Poisson  equation  (18)  is  then  solved  for  the  in¬ 

terior  values  of  x  from  which  (— )  =  and  (— )  =  1^-  are  obtained  ,.t  the 

m  x  9x  m  x  3y 

interior  grid  points.  The  initialized  wind  at  all  interior  grid  points 
is  then  given  by 


V.  =  [(-)  +  <-)  It  +  [(-)  +  (-)  ]} 

in  Lm  ^  m  xj  |_m  1J1  m  xj 


(22) 


(D) .  Boundary  Values 


Data  from  both  the  initial  and  final  time  of  the  forecast  are  pro¬ 
cessed  through  the  preceding  steps.  Then  boundary  values  for  each  prog¬ 
nostic  variable  at  each  location  on  the  exterior  of  the  grid  lattice  are 
calculated  using  a  linear  interpolation  in  time.  These  boundary  values 
are  calculated  for  each  time  step  of  the  forecast.  Time  steps  of  six, 
three,  and  1.5  minutes  are  used  for  the  single,  double,  and  quadruple 
resolution  cases.  In  addition,  the  time  step  increment  (rate  of  change 
of  quantity  at  a  grid  point  multiplied  by  appropriate  time  step  length 
in  seconds)  is  calculated  for  grid  points  located  one  and  two  rows  in  from 
the  boundary  for  each  prognostic  variable  in  order  to  provide  sponge- 
type  boundary  conditions.  The  sponge  boundary  scheme  involves  the  combina¬ 
tion  of  one-third  computed  value  and  two-thirds  increment  value  at  grid 
points  located  one  row  from  the  boundary,  and  two-thirds  computed  value 
and  one-third  tendency  value  at  two  rows  in  from  the  boundary. 
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IV.  Conclusion 


This  report  describes  a  variable  resolution  limited-area  numerical 
weather  prediction  model  developed  to  test  the  effect  on  forecast  accu¬ 
racy  of  changing  spatial  resolution.  Unfortunately,  the  contract  period 
elapsed  before  the  model  could  be  successfully  run  using  atmospheric  data 
Therefore  no  conclusions  can  bo  diawn  concerning  the  relationship  between 
model  resolution  and  forecast  accuracy.  It  is  recommended  that  the  model 
formulated  in  this  effort  be  utilized  to  examine  the  nature  of  this  rela¬ 
tionship  in  finite-difference  models. 


B.  AUTOMATED  WEATHER  SYSTEM 


Under  the  contract  SASC  was  responsible  for  developinq  and  implement¬ 
ing  computerized  data  processing  methods  and  techniques  to  evaluate  new 
or  improved  weather  sensors  and  observational  techniques  and  short-term 
prediction  methods  for  weather  parameters  critical  in  aviation. 

I.  AFGL  Weather  Test  Facility  Support 

SASC  activity  in  this  area  was  largely  based  on  data  produced  by  two 
acquisition  systems  installed  at  the  AFGL  Weather  Test  Facility  (WT1  )  at 
Otis  AFb ,  Cape  Cod.  These  systems  are  the  200-Channel  Doric  Data  L<  gger 
and  the  Modular  Automated  Weather  System  (MAWS).  For  th^m,  SASC  developed 
four  large  CDC  6600  computer  programs:  REQUOTE,  PWP,  MAW.1E,  and  VERIFY. 

In  addition,  many  associated  programs  were  written  to  obtain  specific 
statistical  analyses  and  discrete  data  displays. 

(A) .  DORIC 

The  Doric  Data  Logger  writes  onto  magnetic  tape  voltages  recorded  by 
over  100  weather  sensors  located  on  towers  and  a  ground  site,  const t ucted 
in  a  runway  simulation  configuration. 1  Instruments  included  are  wind  sen¬ 
sors  of  several  types,  temperature  and  dew  point  instruments,  precipita¬ 
tion  gauges,  radiometers,  and  a  variety  of  visibility  meters.  The  boric 
system  gathers  readings  five  times  per  minute,  10  times  in  the  case  of 
wind  sensors,  and  records  the  information  onto  tapes  whi^h  are  shipped 
daily  to  AFGL  at  Hanscom  AFB.  In  the  contract  period  neirly  500  tapes 
were  received  and  processed. 

The  raw  voltage  tapes  are  processed  by  a  computer  program  call'd 
REQUOTE.  This  program  converts  voltages  to  one-minute  mean  meteoro 
logical  unit  data  which  in  turn  are  used  for  display  or  further  analysis. 

Data  received  from  REQUOTE  were  fed  to  the  Present  Weather  Program 
(PWP)  in  order  to  produce  automatic:  "present  weather"  observations  <weu- 

^The  AFGL  Weather  Test  Facility  is  further  described  in  Gerlach,  ed . 

(1980),  pp.  21-21. 
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ther  conditions  and  obstructions  to  vision)  comp.it  at  1  e  to  those  re¬ 
ported  by  a  human  weather  observer.  By  using  a  doe  i  • .  u>n-t  t  ■  e  .ippi.mch 
the  PWP  was  successful  in  discriminating  weather  conditions  such  as  inoder 
ate  or  heavy  drizzle,  rain,  snow,  freezing  drizzle  oi  freezing  rain. 
Obstructions  to  vision  included  blowing  dust  or  sand,  blowing  snow,  smoke 
haze,  fog,  ground  fog,  or  ice  fog. 

Hourly  results  of  automated  observations  weie  c  m;  at  od  by  computer 
plot  with  FAA  human  observations  for  data  collected  over  the  period  March 
1978  -  April  1979.  These  plots  were  used  by  AFGL  in  a  study  on  automa¬ 
tion  of  visual  weather  observations.  The  results  showed  about  84  percent 
agreement.  Lack  of  agreement  for  the  remaining  hours  was  attributed,  in 
part,  to  lack  of  responsiveness  of  the  rain  gauge  to  small  amounts  of  pre 
cipitation  and  in  part  to  the  spatial  separation  of  the  observer  and  the 
test  facility. 

The  FAA  obstruction  to  vision  reports  for  the  same  period,  when  plot 
ted  in  comparison  with  the  results  obtained  from  the  PWP,  showed  similar 
results.  The  difference  between  the  two  methods  was  shown  by  the  AFGL 
study  to  fall  mostly  in  the  five  to  eight  mile  visibility  range.  FAA 
regulations  require  an  obstruction  to  vision  report  only  if  visibility 
is  less  than  seven  miles.  In  the  ranges  close  to  this  limit  variations 
could  occur.  In  general,  however,  analysis  of  the  results  of  the  PWP  us¬ 
ing  the  data  collected  appeared  to  verify  the  approach  taken. 

In  addition  to  the  Present  Weather  Program  model ,  data  were  fed 
from  REQUOTE  to  many  statistical  routines  to  determine  correlations, 
scatter,  and  goodness-of-f it ,  Data  from  REQUOTE  were  also  used  to  evalu¬ 
ate  now  or  improved  weather  sensors.  Wind  sets  such  as  the  Rosemount, 
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J-Tcc,  and  Gill  were  compared  with  the  Climatronics.  Several  models  of 
small-volume  forward-scatter  meters  (Fumosens)  were  evaluated  by  varying 
calibration  equations  and  plottinq  the  conversions,  in  extinction  coeffi¬ 
cients,  in  comparison  with  transmissometer  data. 

Another  imixartant  use  of  the  data  base  was  in  the  slant  visual  range 
prediction  model  discussed  below  (pp.  27,  30).  To  determine  probability 
estimates  for  below-limits  SVR  conditions,  SASC  wrote  routines  to  incor¬ 
porate  Markov  and  Equivalent  Markov  techniques  into  REQUOTE.  Also,  a 
plot  program  was  written  to  display  the  predicted  probability  of  below- 
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See  fieri ach,  ed .  (1980),  pp.  20-29,  for  a  detailed  description  of  the 

wind  set  tests. 
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limit  visibility  using  the  Markov  and  Equivalent  Markov  techniques. 

A  Doric  wind  shear  routine,  DORWIND,  was  developed  1  y  SASC  to  evalu¬ 
ate  the  spatial  homogeneity  of  wind  measurements  collected  at  towers  A,  y>, 
and  X  at  WTF.  Sensors  located  at  various  tower  levels  rijjort  values  used 
to  compare  wind  shears  and  statistics  such  as  mean,  standard  deviation, 
standard  error,  root  mean  square  error,  bias,  correlation,  and  line  of 
best  fit.  DORWIND  is  run  on  a  regular  basis  with  each  Doric  tape. 

(B) .  MAWS  (Modular  Automated  Weather  System) 

A  second  data  acquisition  system,  MAWS,  was  installed  at  the  WTF  in 
late  Summer,  1979.  Using  the  same  instruments  as  the  Do)ic,  with  the  addi¬ 
tion  of  three  rotating  beam  ceilometers  for  determining  <  loud  base  height, 
the  MAWS  records  one-minute  means  in  meteorological  unit:  onto  magnetic 
tape.  It  also  has  the  capability  of  interspersing  raw  voltages  with  the 
weather  data  means  to  provide  test  information.  In  addition,  three  remote 
sites  on  Cape  Cod,  instrumented  to  give  additional  information  particu¬ 
larly  on  sea  fog,  are  tied  into  the  MAWS  retrieval  system.  The  sensor  in¬ 
formation  contained  in  the  data  records  is  relayed  over  telephone  long 
lines  to  Hanscom  AFB  and  packed  onto  tapes  at  AFGL.  Because  of  packinq 
one-minute  mean  data,  rather  than  12-second  voltaqe  data  as  in  Doric, 

MAWS  tapes  hold  15  to  20  days  of  data  per  tape. 

Computer  proctram  MAWSE  was  developed  by  SASC  to  accomplish  the  un¬ 
packing,  processing,  and  display  of  data  acquired  from  the  sensors  at 
Otis.  An  editing  routine  was  devised  to  recover  data  records  sometimes 
marked  by  faulty  identifications  and  other  errors  transient  in  the  system, 
and  was  implemented  in  MAWSE. 

Using  MAWS  data,  program  VERIFY  was  developed  by  SASC  to  perform  a 
thorough  verification  analysis  of  short-range  probability  and  categorical 
forecasts.  Probability  forecasts  are  substantiated  using  rank  probability 
scores.  Brier  scores,  a  "T"  statistic,  a  measure  for  skill,  and  a  Fried¬ 
man's  Analysis  of  Variance  test.  Categorical  forecasts  are  verified  us¬ 
ing  a  Chi-Square  statistic,  percent  correct,  threat  score,  and  bias. 

VERIFY  also  updates  a  master  file  containing  raw  statistics  from  all  epi¬ 
sodes  of  low  ceiling. 

Two  routines,  WINDCOM  and  VSHYCOM,  Wert;  written  by  .'•ASC  to  compare 


statistically  wind  and  visibility  observations  rcjjorted  by  MAWS  (at  Scott 
AFB)  and  a  human  observer  at  Scott  durinq  the  period  Jan.  1977  -  Jan.  197 3. 
MAWS  values  were  averaged  over  a  ten-minute  period  and  compared  with  a  human 
observation  taken  once  durinq  the  same  ten-minute  interval .  WINDCOM  and 
VS BY  'OM  computt  means,  standard  deviations,  variances,  bias,  correlation, 
root  mean  squaie  error,  percent  root  mean  square  error,  standard  error, 
covariance,  anc.  lino  of  best  fit.  Statistics  were  updated  as  successive 
tapes  were  processed. 

(C) .  Other 

Other  software  written  by  SASC  included  a  data  retrieval  program  to 
process  tapes  from  the  Monitor  Labs  Data  Acquisition  System  used  for  test¬ 
ing  laboratory  instruments.  This  program  provides  interpretation  of  any 
number  of  data  channels  up  to  200,  and  any  number  of  scans  of  the  channels 
capable  of  fitting  on  a  single  tape.  Also  devised  was  a  data  reduction 
program  for  the  Particle  Measuring  System  Inc.  spectrometer,  which  measures 
fog  droplet  size. 

II.  Development  of  Observational  Techniques 

(A) .  AWDS  Scenario  Software 

To  meet  requirements  for  accurate  and  timely  reporting  of  airfield 
weather  conditions,  the  Air  Weather  Service  (AWS)  plans  to  modernize  its 
Base  level  weather  support  in  the  1980's.  The  Automated  Weather  Distribu¬ 
tion  System  (AWDS)  will  encompass  the  modernization  of  Air  Force  weather 
station  functions.  AWDS  will  serve  as  a  data  products  handler,  producing 
computer-generated  alphanumeric  and  videographic  displays  of  meteorologi¬ 
cal  information.  Critical  weather  parameters  which  affect  airfield  opera¬ 
tions  include  visibility,  wind  speed  and  direction,  turbulence,  and  precipi¬ 
tation.  AWDS  will  provide  the  technology  with  which  these  parameters  will 
be  collected,  analyzed,  and  displayed. 

SASC  meteorologists  developed  software  applicable  to  the  AWDS  sce¬ 
nario  through  adaptation  and  exploitation  of  software  developed  for  the 
Man-computer  Interactive  Data  Access  System  (McIDAS) .  Development  concen¬ 
trated  on  routines  to  display  standard  and  derived  meteorological  para¬ 
meters.  These  new  routines  access  current  McIDAS  library  moiules  to  pro- 
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duce  both  plotted  and  analyzed  weather  data. 

Specific  operational  sottware  includes  a  proqram  which  plots  upper 
air  wind  observations  on  constant  height  levels.  This  program,  WIH1IGT, 
displays  either  wind  vectors,  wind  flags,  or  numerical  values  of  wind 
speed  and  direction.  A  total  of  16  levels,  from  1,000  to  50,000  feet,  may 
be  plotted  on  base  maps  ranging  from  full  USA  to  individual  states.  De¬ 
rived  atmospheric  stability  parameters  and  estimates  of  precipitable  water 
are  also  available.  This  program,  STABLE,  provides  the  option  of  display¬ 
ing  either  plotted  or  analyzed  information.  Quantities  available  for  dis¬ 
play  include  K-index,  totals  index,  lifted  index,  convective  temperature, 
and  precipitable  water.  Information  may  be  displayed  on  Mercator  projec¬ 
tions  or  overlaid  on  satellite  imagery.  Output  can  be  directed  to  either 
a  CRT  or  line  printer. 

(B) .  Cloud  Base  Height  Estimation 

Data  storage  and  manipulation  routines  were  developed  to  perform  the 
initial  evaluation  of  a  series  of  algorithms  used  to  translate  rotating 
beam  ceilometer  (P.BC)  measurements  into  estimates  of  cloud  base  height 
(from  0  to  6,000  ft)  and  total  sky  condition  (clear,  scattered,  broken,  or 
overcast) .  A  total  of  15  different  algorithms  was  tested. 

Preliminary  results  suggested  modifications  to  the  cloud  base  height/ 
sky  condition  program  CLOUD.  Originally,  CLOUD  produced  output  which  in¬ 
cluded  sky  condition  estimates  for  up  to  five  different  layers  between  the 
surface  and  8,700  ft.  RBC  measurement  error  (inherent  in  the  instrument) 
necessitated  lowering  the  height  of  the  maximum  cloud  layer  from  8,700  to 
6,000  ft.  With  this  drop  in  maximum  cloud  height  the  maximum  number  of 
cloud  layers  was  also  reduced  (from  five  to  three)  to  produce  a  more  real¬ 
istic  depiction  of  the  sky  condition. 

In  all  a  total  of  14  episodes  (62  individual  data  sets)  was  processed 
during  the  period  from  June  through  September,  1980. 

Figure  B- 1  is  a  flow  chart  of  the  storaqe  and  manipulation  software. 
Each  step  represents  a  program  used  either  to  reduce  or  store  processed 
information.  The  individual  programs  are  as  follows: 

TRANS  -  Input  to  this  program  comes  directly  from  CLOUD.  The  informa¬ 
tion  is  coded  and  packed  to  speed  processing  and  reduce  stor- 
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Figure  B-l.  Data  storage  and  manipulation  software. 
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age.  Output  is  stored  in  separate  episodes  containing  as 
many  as  ten  different  methods  (calculations)  for  evaluation. 

STORE  -  Used  to  store  the  packed  data  from  TRANS  on  magnetic  tape. 
Information  is  filed  in  separate  episodes. 

HGTPLT  -  Produces  time  vs.  height  ink  plots  of  any  one  method  (1  - 
15)  for  a  particular  episode.  Plots  represent  altitude  of 
individual  cloud  decks,  the  condition  of  each  layer,  whether 
the  layer  is  continuous  with  time;,  the  heiqht  of  measured 
ceiling,  and  the  presence  of  a  layer  above  the  measured  ceil¬ 
ing.  Figure  B-2  is  a  sample  of  HGTPLT  output. 

FAADTA  -  Produces  C'DC  6600  data  files  from  sky  condition  reports 
taken  by  FAA  observers.  The  information  is  coded  and  key 
punched  in  a  form  compatible  with  the  information  received 
from  the  automatic  routine  and  provides  ground  truth  for  the 
evaluation. 

SORT  -  Produces  a  single  data  file  by  combining  both  the  observed  data 
set  and  calculated  estimates. 

COMPAR  -  Produces  computer  ink  plots  of  observed/ca Lculated  ceiling 
heights  on  time  vs.  height  axes  for  any  one  method  during  a 
particular  episode.  Figure  B-3  is  a  sample  of  COMPAR  output. 

DTADMP  -  Calculates  statistical  information  used  to  evaluate  the 

cloud  height  and  sky  conditions  algorithms.  For  every  method 
tested,  DTADMP  tabulates  observed  vs.  calculated  number  of 
layers,  total  sky  cover,  and  ceiling  reports,  and  calculates 
the  following  statistics: 

Percent  Correct 

Percent  Agreement 

Probability  of  Detection  (POD) 

Threat  Score  (TS) 

False  Alarm  Ratio  (FAR). 

Method  statistics  are  presented  for  all  past,  present,  and 
total  (past  +  present)  episodes. 

III.  Development  of  Short-Range  Prediction  Techniques 

(A).  Objective  Prediction  of  Slant  Visual  Range 

BASC  meteorologists  aided  AFG]  scientists  in  developing  an  objective 
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technique  for  specifying  visibility  conditions  below  I-'AA  Category  1  and 
II  aircraft  landing  guidelines  from  forward  scatter  meter  (FSM)  measure¬ 
ments.  Atmospheric  extinction  coefficients  measured  by  an  array  of  FSM's 
on  three  instrumented  towers  at  AFGL's  Weather  Test  Facility  were  used  to 
specify  i  he  Slant  Visual  Range  (SVR)  in  the  landing  approach  zone. 

The  accuracy  and  reliability  of  three  prediction  techniques  used  to 
estimate  SVR  a  .  forecast  intervals  of  2,  5,  10,  30,  and  6C  minute;:,  were 
tested  aid  evaluated.  Techniques  consisted  of  an  adaptation  of  a  Markov 
chain  knawn  as  the  "Ornstein-Uhlenbeck"  process.  Regression  Estimation  of 
Event  Pr labilities  (REEP) ,  and  an  Equivalent  Markov  method  based  on  the 
classical  Markov  transition  matrix.  Prediction  equations  were  formulated 
using  twenty  cases  of  reduced  visibility  (6,300  minutes  of  data).  An 
independent  data  set  of  twenty-five  episodes  (9,900  minutes  of  data)  was 
used  to  compare  the  techniques. 

Results  of  this  study  indicated  that  in  cases  of  reduced  visibility 
caused  by  advection  fog,  a  prediction  scheme  which  utilizes  SVR  measure¬ 
ments  yields  far  more  accurate  forecasts  of  below-limit  conditions  than 
the  stan lard  Runway  Visual  Range  (RVR)  measurements  being  used.  Proba¬ 
bility  of  detecting  reduced  visibility  conditions  for  both  Category  I  and 
II  landing  operations  is  91  percent  or  greater. 

Of  the  three  prediction  techniques  evaluated,  the  Equivalent  Markov 
scheme  c demonstrated  a  higher  degree  of  reliability  and  accuracy  than  did 
the  Mar) ov  or  REEP  techniques. 

A  paper  entitled  "Objective  Prediction  of  Slant  Visual  Range  During 
Advective  Fog  Situations,"  which  documents  the  results  of  this  study  (co¬ 
authored  by  SASC  and  AFGL  meteorologists) ,  was  presented  at  the  American 
Meteorological  Society's  Eighth  Conference  on  Weather  Forecasting  and 
Analysis,  Denver,  June  10  -  13,  1980.  The  paper  appears  in  the  preprints 
of  the  conference.^ 

(B) .  Objective  Prediction  of  Wind 

Study  of  mesoscale  wind  variability,  leading  to  the  development  of 
automated  wind  analysis  and  prediction  methods,  required  establishment  of 

Geisler,  E.B. ,  and  R.  C.  Dengel,  1980:  Objective  prediction  of  slant 
visual  ranqe  d zring  advective  fog  situations.  Preprint,  Eighth  Conference 
on  Weather  For ^casting  and  Analysis,  AMS,  Boston;  pp.  420-423. 
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a  suitable  data  base.  Data  tapes  generated  by  the  AFGL  Mesonetwork  were 
selected  for  this  purpose.  The  Mesonet  collected  meteorological  data 
(wind,  temperature,  dew  point,  visibility,  and  rainfall)  from  September 
1972  through  May  1976  at  26  specially  designed  weather  stations  in  north¬ 
eastern  Massachusetts . 

A  three-phase  project  was  established  to  edit  the  data,  generate  a 
climatology,  and  perform  statistical  analyses  aimed  at  developing  an  ob¬ 
jective  wind  prediction  technique. 

Programs  designed  by  SASC  to  edit,  condense,  and  group  these  data  into 
a  more  complete  data  set  are  MESOFV,  HOURAVE,  HSTOGRM ,  and  EDIT. 

MESOPV  converts  voltage  values,  taken  at  ten-second  intervals  from 
each  of  the  26  Mesonet  stations,  into  one-minute  meteorological  means. 

These  one-minute  values  of  wind  direction  and  speed,  temperature,  dew 
point,  and  extinction  coefficient  (exco)  are  written  onto  a  permanent 
file.  This  file  eliminates  further  processing  of  the  Mesonet  magnetic 
tape,  thus  decreasing  computer  processing  time. 

One-minute  means  are  then  computed  into  hourly  and  daily  averages 
using  the  program  HOURAVE.  Hard  copies  of  these  averages,  along  with 
maximums,  minimums,  and  standard  deviations,  are  printed  for  each  Mesonet 
file  generated  by  MESOFV. 

HSTOGRM  was  designed  as  a  supplement  to  HOURAVE,  to  make  easier  the 
identification  of  spurious  data.  This  program  generates  a  histogram  of 
each  variable  at  each  station  for  the  period  covered  by  one  Mesonet  tape. 

In  addition,  a  list  of  variables  and  times  in  which  an  erroneous  value 
seems  to  occur  is  printed  at  the. end  of  each  program. 

MESOFV  is  run  only  once  with  HOURAVE  and  HSTOGRM  for  initial  data 
processing  purposes.  All  other  editing  is  done  using  one-minute  means 
from  the  permanent  file  generated  by  MESOFV,  hereafter  referred  to  as 
netf ile. 

Hard  copies  produced  by  HOURAVE  and  HSTOGRM  are  analyzed  and  spurious 
data  identified.  A  set  of  data  cards  is  punched  containing  times  and  sta¬ 
tions  at  which  a  variable  has  been  flagged  as  bad.  These  cards  are  used 
as  input  for  EDIT,  which  searches  netfile  for  flagqed  variables  and  re¬ 
places  them  with  a  value  of  999.  Whenever  this  value  is  encountered  by 
later  programs,  it  is  ignored.  EDIT  also  insures  that  there  are  no  time 
gaps  on  the  file,  by  filling  in  missing  minutes  and  giving  the  correspond¬ 
ing  variables  a  value  of  999.  When  edging  is  complete,  netfile  is  re- 


placed  by  an  edited  version  of  the  one-minute  means.  HOURAVE  and  HSTO- 
GRM  then  produce  updated  versions  of  their  respective  outputs  to  be  re¬ 
analyzed  for  further  editing. 

This  editing  process  is  repeated  (usually  three  or  four  times)  until 
a  fully  edited  netfile  exists.  The  edited  netfile  is  later  packed  by  pro¬ 
gram  MAGFILE  onto  a  new  magnetic  tape  (Mesoxx) .  Ten  netfiles  are  placed 
on  one  Mosoxx  tape.  Eventually  these  one-minute  mean  packed  tapes  will 
replace  the  Mesonet  tapes. 

Each  of  these  programs  was  designed  to  make  intensive  data  reduction 
as  selective  as  possible.  Because  of  the  amount  of  processing  time  and 
data  involved,  job  turnaround  time  was  often  slow.  To  contract  completion, 
40  Mesonet  tapes  were  fully  edited  and  four  new  magnetic  tapes  created 
(MesoOl  -  Meso04) .  With  this  newly  created  data  base,  preparations  for 
analysis  of  mesoscale  wind  variability  were  begun. 

SAS<’  prepared  a  series  of  programs  to  generate  a  climatological  atlas 
of  monthly,  seasonal,  and  annual  weather  at  various  times  of  the  day.  One 
of  these  programs  is  MAINFRE,  which  categorizes  wind  direction  and  speed, 
temperature,  dew  point,  and  exco  values,  and  stores  them  on  permanent  ran¬ 
dom  access  files.  These  files  are  used  to  generate  tables  containing  rela¬ 
tive  frequency  distributions,  maximum  and  minimum  values,  and  means  of 
given  periods  of  time.  SASC  also  wrote  programs  to  print  these  tables  on 
hard  copy. 

The  tab.es  contained  in  the  atlas  are: 

Wind  direction  vs.  wind  speed 

Wind  direction  vs.  exco 

Wind  direction  vs.  temperature 

Wind  speed  vs.  temperature 

Wind  speed  vs.  exco 

Maximum  temperatures 

Minirum  temperatures 

Mean  temperatures 

Maxirum  wind  speed  and  direction 

Exco 

Each  frequency  distribution  is  dependent  on  station,  month,  and  time 
factors.  This  climatological  .summary  is  run  and  updated  with  each  Mesoxx 
tape.  When  th<  atlas  is  completed,  results  will  be  interpreted  for 
effects  of  top<  qraphy,  diurnal  change,  and  synoptic  scale  weather  patterns. 


C.  SATELLITE  METEOROLOGY:  OBJECTIVE  PROCESSING  OE  SATELLITE  IMAGEI Y 


I.  Identification  of  Cloud  Types 

SASC  scientists  investigated  an  objective  satellite  imagery  classi¬ 
fication  technique  which  identifies  the  cloud  type  of  a  visible  or  infra¬ 
red  image  sample  25  x  25  n  mi  in  size.  Defense  Meteorological  Sate! life 
Program  (DMSP)  2/3  n  mi  visible  and  infrared  data  from  an  area  in  the  vi¬ 
cinity  of  the  Mediterranean  Sea  were  used  for  this  study.  These  data  on 
tape  were  supplied  to  AFGL  by  the  Air  Force  Global  Weather  Central  (AFGWC) 
The  classifier,  a  modification  of  an  earlier  technique  applied  to  visible 
data  only,  extracts  information  about  the  texture  and  patterns  of  clouds 
in  an  imagery  sample  by  examining  the  coefficients  of  its  discrete  Four¬ 
ier  transform.  It  is  capable  of  classifying  a  37  x  37-pixel  (25  x  35  n 
mi)  image  in  less  than  one  second  on  AFGL’s  CDC  6600  computer.  The  classi 
fier  was  tested  on  a  total  of  143  image  samples,  and  was  allowed  to  choose 
the  cloud  type  of  an  image  sample  from  one  of  nine  types:  stratocumulus , 
stratus,  cumulus,  clear,  altocumulus,  altostratus/nimbostratus ,  cirrus, 
cirrostratus ,  and  cumulonimbus.  On  average,  7  to  8  image  samples  were 
correctly  identified  out  of  every  10  it  was  allowed  to  classify.  A  scien¬ 
tific  report,  "Performance  of  the  Discrete  Fourier  Transform  Satellite 
Imagery  Classification  Technique,"  submitted  for  publication  in  June  1980, 
offers  a  complete  description  of  this  study.1 

II.  Snow/Cloud  Discriminator 

In  June  1979  the  DMSP  launched  an  experimental  snow/cloud  disc  rimina- 
tor  sensor  (SSC)  in  a  polar  orbit.  The  sensor  is  a  linear  array  ot  detec¬ 
tors  sensitive  to  reflected  sunliqht  in  the  1.51  -  1.63  micron  spectral 
region.  Visible  (0.4  -  1.1  micron),  infrared  (10  -  12  micron) ,  anc  snow/ 
cloud  data  for  December  1979  were  sent  on  tape  to  AFGL  by  AFGWC.  <  overage 
was  limited  to  near  local  noontime  orbits  over  eastern  North  and  C’t  ntral 
America,  and  extreme  northern  South  America.  Because  ot  (lower  supply  prob 

^d ' Entremont ,  R.P.,  1980:  Performance  of  the  discrete  Fourier  transform 
satellite  imagery  classification.  AFGt-tr-80-0 1 75 .  Scientific  Refiort 
No.  3,  Contract  F  19628-79-C-003  1 . 
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lems,  no  data  wore  available  alter  December  1979. 

Snow/cloud  imaqery  contains  useful  information  with  which  cloud/no 
cloud  and  snow/  10  snow  decisions  can  be  more  confidently  made,  inasmuch  as 
the  snow/cloud  ipectral  reqion  exhibits  several  unique  physical  char.ictcr  - 
istics.  In  thi ?  reqion  of  the  near  1R,  water  droplets  are  good  reflector  s 
of  sunliqht,  ic  i  particles  are  not  so  qood,  and  oceans,  lakes,  ana  ot..er 
bodies  of  water,  alonq  with  stv ..w-co'. or  'd  ground,  .ire  poor  :e:io :  r.  . 

When  land  is  not  covered  by  snow  or  ion,  it  appears  briqhtor  in  t no  ;  ,aw, 
clou  1  bandwidth  than  it  does  in  the  visible.  Coastlines  and  shores  thus 
contrast  well  with  surrounding  waters,  which  appear  black. 

In  the  visible  and  in  the  10  -  12  micron  IR,  snow,  low  stratiform, 
and  low  cumuliform  clouds  tend  to  appear  brighter  than  clear  land  or  ocean. 
Thus  fog  and  low  stratus  can  be  confused  with  snow-covered  ground.  How¬ 
ever,  snow-covered  land  is  dark  in  the  snow/cloud,  while  low  stratiform 
and  cumuli  form  clouds  composed  of  water  droplets  are  bright.  The  sharp 
distinction  between  snow  (dark)  and  liquid  clouds  (bright)  earns  the  SSC 
the  name  "snow/cloud  discriminator." 

Snow/cloud  data  also  can  help  satellite  imagery  analysts  more  confi¬ 
dently  determine  whether  clouds  are  composed  of  ice  particles  or  water 
droplets.  Clouds  that  are  bright  in  the  visible  and  warm  in  the  IR  appear 
bright,  in  the  snow/cloud  (water  droplets)  ,  while  clouds  that  are  bright  in 
the  visible  but  cold  in  the  IR  appear  gray  in  the  snow/cloud  (ice  particles) . 

Figure  C- 1  shows  the  geographical  locations  of  the  images  in  Figures 
C-2  and  C-3.  I  hose  figures  contain  SSC  imagery  which  illustrates  some  of 
the  features  discussed  in  the  preceeding  paragraphs.  Figure  C-2  contains 
visible,  IR,  and  snow/cloud  imagery  over  a  region  of  northwestern  South 
Ame. ica.  In  the  lower  part  of  each  image  a  large  group  of  cumulonimbus 
cloi ds  can  be  seen.  These  clouds,  between  250  and  300  n  mi  wide,  appear 
bricht  in  the  visible  and  cold  in  the  IR.  Snow/cloud  imagery  in  the  corres- 
ponoing  bright/ cold  areas  of  the  visible/IR  pictures  is  gray,  indicating 
that  the  cloud  tops  consist  of  ice  particles.  Surrounding  lower-level 
waf'  r  clouds  (extreme  lower  right)  are  much  brighter. 

Figure  C-3  is  imagery  over  Cape  Cod,  the  coast  of  Maine,  and  Nova 
Sco  ia.  The  blight  clouds  (warm  in  the  IR)  in  the  lower  right  of  each 
vis  ble  image  ere  st ratocumulus .  The  fact  that  these  clouds  appear  bright 
in  he  snow/cloud  indicates  that  they  consist  of  water  droplets.  Note 
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that  the  Maine  and  Nova  Scotia  coastlines  stand  out  against  the  ocean  in 
the  snow/cloud  better  than  they  do  in  the  visible.  Unfortunately,  snow 
cover  in  eastern  North  America  was  very  sparse  in  December  1979  when  these 
data  were  collected.  Thus  there  were  not  many  good  snow-covered  ground 
imagery  samples  to  use  as  illustrations. 

Each  image  in  Figures  C-2  and  C-3  is  labeled  with  its  sensor  and  reso¬ 
lution  at  satellite  subpoint.  The  cross-track  resolution  of  a  visible  or 
infrared  pixel  is  1.66  n  mi  at  satellite  subpoint,  and  decreases  to  2.95 
n  mi  at  scan  edge.  Along-track  resolution  is  close  to  3.0  n  mi.  The 
visible  and  IR  images  to  either  side  of  the  6  n  mi  snow/cloud  image  were 
generated  by  averaging  the  3  n  mi  visible  and  IR  data  in  order  to  make 
them  look  like  (match  up  with)  the  snow/cloud  imagery.  (This  was  done  by 
AFGWC,  but  a  similar  set  of  averaged  visible  and  IR  images  was  also  gen¬ 
erated  by  SASC  at  AFGL . )  The  cross-track  resolution  of  the  coarser  snow/ 
cloud  data  is  4.1  n  mi  at  satellite  subpoint,  decreasing  to  8.3  n  mi  at 
scan  edge.  Along-track  resolution  of  SSC  decreases  from  6.6  n  mi  at  sub- 
point  to  9.1  n  mi  at  scan  edge. 

When  visible  and  IR  images  of  3  n  mi  resolution  are  averaged 
(smoothed)  to  correspond  pixel  for  pixel  to  snow/cloud  imagery,  small- 
scale  variabilities  in  cloud  brightness  and  temperature  are  lost.  As  a 
consequence,  a  significant  amount  of  small-scale  cloud  feature  informa¬ 
tion  is  lost.  Knowing  the  "footprint"  sizes  (resolutions  of  each  pixel 
edge)  and  geographical  locations  of  each  3  n  mi  visible  and  IR  pixel  and 
each  6  n  mi  snow/cloud  pixel  enables  the  higher  resolution  visible/IR 
data  to  be  mapped  into  the  coarser  snow/cloud  imagery.  Thus  the  several 
visible  and  IR  3  n  mi  pixels  that  lie  within  the  same  area  as  a  6  n  mi 
snow/cloud  pixel  can  be  found.  More  information  about  the  type  of  feature 
an  SSC  pixel  represents  can  be  obtained  by  observing  the  smaller-scale 
brightness  variations  of  those  3  n  mi  pixels.  Such  analyses  can  help  an 
automated  classifier  determine  the  cloud  or  surface  type  of  an  image 
sample. 

Figure  C-4  contains  graphical  displays  of  the  relationships  among 
averaged  visible,  averaged  IR,  and  snow/cloud  grayshades.  These  graphs 
were  constructed  by  randomly  choosing  an  interesting  cloud  feature  on  an 
image,  subjectively  classifying  it,  and  then  using  a  McIDAS  command  to 
extract  the  visible,  IR,  and  SSC  graysnade  values  of  that  feature.  Re- 
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suits  were  manually  tabulated.  These  graphs  point  out  grayshade  thres¬ 
holds  among  different  meteorological  situations. 

In  Figures  C'-4a  and  C-4b  it  is  seen  that  the  range  between  land  and 
ocean  backgrounds  in  the  snow/cloud  is  greater  than  it  is  in  either  the 
visible  or  the  IR.  Note  how  well  water  clouds  are  separated  from  sr, iw- 
covered  ground  in  the  SSC,  and  also  how  well  water  clouds  are  separated 
from  ice  clouds;.  The  thresholds  illustrated  in  these  graphs  could  be 
used  as.  input  to  a  cloud  type-  and  phase-determining  algorithm. 

Beyond  the  curves  bounding  the  cloud  types  in  the  graphs  of  Figure 
C-4 ,  no  grayshades  were  found  in  the  images  that  were  examined.  It  is 
interesting  to  note  in  Figures  C-4a  and  C-4b  that  the  SSC  is  often  satur¬ 
ated,  primarily  by  thick,  dense  water  droplet  clouds.  In  the  visible  and 
IR,  sensor  saturation  would  most  likely  be  caused  by  thick  cirrostratus 
or  well  developed  cumulonimbus. 

In  summary,  certain  characteristics  of  snow/cloud  imagery  logically 
determine  how  best  to  use  it  in  an  automated  mode.  Sharp  snow/bare 
ground  and  ice  cloud/water  cloud  contrasts  dictate  that  an  automated  snow/ 
cloud  imagery  processing  algorithm  should  distinguish  snow  on  the  ground 
from  bare  ground  or  low  clouds  and  identify  cloud  phase.  Other  potential 
applications  of  those  data  include  determination  of  cloud  type,  frac¬ 
tional  cloud  cover,  cloud  top  height,  thickness,  and  temperature,  and 
estimations  of  mesoscale  visibility. 

III.  Brightness  Variation  Analysis 

In  an  effort  to  utilize  geostationary  satellite  data  for  determin¬ 
ing  cloud  cover  amounts,  SASC  scientists  developed  and  implemented  soft- 

2 

ware  to  analyze  ODIS  (Offline  Data  Ingest  System)  data.  ODIS  data  con¬ 

sist  of  one-mile  resolution  hourly  images  centered  on  the  northeastern 
US  as  received  from  the  GOES  and  SMS  series  of  geostationary  satellites. 
Data  from  the  period  March  1977  -  April  1979  were  processed.  Seven  sites  - 
Syracuse,  Buffalo,  Wilkes-Barre,  Philadelphia,  Atlantic  City,  Lynchburg, 
and  Wallops  Island  -  were  selected  as  study  stations  to  provide  a  geo¬ 
graphical  mixture  of  the  samplinq  region.  Hourly  weather  observations 

2 

See  Gerlach,  ed.  (1980)  for  detailed  description  of  the  data  reduction 
method  and  software;  pp.  74  -  81. 
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from  these  stations  were  used  as  the  ground  truth  in  the  study. 

The  first  step  in  the  analysis  was  a  data  reduction  to  provide  more 
easily  studied  parameters.  To  this  end  program  CENTER  was  developed. 

CENTER  identifies  and  corrects  garbled  or  missing  scan  lines  of  data  and 
bad  row  count  numbers.  Provisions  were  also  made  to  ignore  occasional 
half-hourly  data  which  have  no  corresponding  ground  truth  data.  Then  us¬ 
ing  basic  and  fine-tuned  navigation  information  provided  by  the  McIDAS, 
CENTER  isolates  a  27  x  27  picture  element  array  centered  on  a  particular 
station.  This  array  is  further  divided  into  nine  9  x  9  arrays  and  the 
brightness  count  for  each  picture  element  is  converted  to  a  reflectivity 
value  normalized  for  solar  zenith  angle.  The  average  reflectivity,  maxi¬ 
mum  and  minimum  values,  and  the  standard  deviation  and  range  of  values 
for  each  9x9  array  and  the  entire  27  x  27  array  are  then  calculated  and 
saved.  These  data  must  then  be  checked  to  insure  that  the  navigation  para¬ 
meters  are  applied  to  the  proper  image.  Bad  GDIS  data  headers,  multiple 
headers,  and  incorrect  recorder  start-up  times  occurred  and  the  necessary 
modifications  were  made  to  correct  the  input  data  for  those  tapes. 

The  collected  data  were  then  stratified  into  eleven  categories  of 
opaque  cloud  cover:  clear,  10  percent  opaque  cloud  cover,  20  percent,  30 
percent,  etc.,  usinq  the  ground  truth  data  as  the  basis  for  stratification. 
Program  ALBZGT  performed  a  first  order  statistical  analysis  of  these  data. 
Erroneous  samples  were  checked  for  poor  navigation  and  corrected  or  re¬ 
jected  depending  on  individual  circumstances.  The  eleven  categories  were 
further  subdivided  into  two  classes,  water  clouds  and  ice  clouds,  and  then 
re-analyzed  by  ALBZGT.  The  results  generated  could  be  used  to  determine 
the  characteristics  which  define  each  cloud  cover  class. 

In  a  related  effort  program  THRESH  was  developed  by  SASC  to  aid  in 
constructing  a  cloud  cover  truth  table.  Using  the  same  27  x  27  array  of 
reflectivity  values  isolated  by  CENTER,  program  THRESH  calculates  the  num¬ 
ber  of  picture  elements  in  the  array  that  are  above  a  given  threshold 
value.  Also  calculated  is  the  percent  of  the  total  array  that  number 
represents.  Several  threshold  values  can  be  input  to  determine  which  one 
provides  the  best  ground/cloud  discriminator  value.  Once  this  best  thres¬ 
hold  value  is  determined,  the  percent  of  the  27  x  27  array  containing  clouds 
can  be  correlated  to  a  percent  cloud  cover  as  determined  by  the  ground 
truth  data  producing  the  cloud  cover  truth  table. 


41 


Software  was  also  developed  by  KAKC  so  i  out  i  sts  to  analyze  the  intru- 
red  data  which  correspond  to  the  visual  data  that  were  already  processed . 
Program  ANALYZE  performs  the  same  functions  as  CENTER  except  for  use  of 
thermal  counts,  not  reflectivity  values.  However,  since  the  JR  data  have 
a  resolution  of  two  miles  in  the  east-west  direction  and  four  miles  in 
the  north-south  direction,  the  additional  steps  of  dividing  each  JR  pic¬ 
ture  element  into  eight  one-mile  picture  elements  and  tiier.  aligning  these 
expanded  data  WLth  the  previously  analyzed  visual  data  must  be  performed. 
Unfortunately,  the  visual-IR  alignment  changes  with  each  satellite,  creat¬ 
ing  a  problem  which  remains  to  be  solved. 

IV.  Navigation  Fine-Tuning 

The  earth  location  of  points  on  satellite  images  is  characteristically 
in  error  by  several  kilometers  because  of  inherent  inaccuracies  in  the  nav¬ 
igation  procedure  employed  by  McIDAS.  While  errors  of  this  magnitude  are 
tolerable  for  many  qualitative  purposes,  research  applications  frequently 
require  greater  precision.  A  fine-tuning  method  was  developed  by  SASC  to 
determine  correction  factors  which  are  applied  to  navigated  data  to  achieve 
the  accuracy  required. 

Visible  channel  imagery  recorded  by  the  ODIS  is  read  into  McIDAS  and 
displayed  on  a  TV  screen.  The  positions  of  one  to  four  landmarks  (earth 
features  whose  locations  are  known  to  an  accuracy  greater  than  or  equal 
to  the  resolution  of  the  satellite  sensor)  are  measured.  The  differences 
between  the  locations  of  the  landmarks  as  predicted  by  the  original  navi¬ 
gation  and  their  actual  measured  positions  are  used  to  determine  two  off¬ 
sets,  one  for  sun  line  number  and  one  for  element  number  (roughly  corres¬ 
ponding  to  nort-i/south  and  east/west  offsets) .  These  values  are  used  to 
compute  correction  factors  for  the  navigation.  This  technique,  when  ap¬ 
plied  to  each  inage  time,  is  capable  of  reducing  earth  location  errors  to 
one  kilometer  or  less. 
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D.  Me IDAS  OPERATION  AND  ENHANCEMENT 


I.  Hardware  Acquisitions 

McIDAS  hard  copy  capability  was  improved  by  acquisition  of  a  Model 
146  Memory  Interface  manufactured  by  Leonessa  Engineering  Company  of 
Bedford,  MA.  This  memory  inter face  stores  a  satellite  image  600  lines 
by  672  elements  or  a  total  of  316,000  pixel  elements.  The  interface  ac¬ 
cepts  data  at  a  high  transmission  rate  from  the  computer,  stores  them  in 
memory,  then  re-transmits  the  data  at  a  slower  rate  required  by  the 
Laserfax  Model  850  to  produce  the  final  copy.  The  resulting  8*5"  x  7*j" 
print  is  of  markedly  higher  quality  than  that  produced  by  the  previously 
used  hard  copy  system.^  The  Leonessa  system  also  features  a  grayscale 
adjustment  panel  which  permits  various  enhancements  to  be  set  into  an 
image . 

During  1980  a  Model  147  Radar  Interface,  also  from  Leonessa 
Engineering  Company,  was  purchased  to  add  to  the  existing  McIDAS  com¬ 
ponents.  A  memory  telephone  with  19  pre-selectable  numbers,  each  corres¬ 
ponding  to  a  particular  radar  site,  was  incorporated  with  the  Model  147 
into  the  McIDAS  Facsimile  Interface.  The  resulting  system  allows  the 
user  to  call  up  and  display  in  real-time,  radar  data  for  sites  located 
in  the  US.  The  19  pre-selected  radar  sites  may  be  changed  according  to 
the  user's  requirements. 

An  engineering  study  was  initiated  to  determine  what  components 
are  needed  to  enable  McIDAS  to  receive  GOES  satellite  data  from  the  experi¬ 
mental  VISSR  Atmospheric  Sounder  (VAS)  satellite  launched  in  September 
1980.  The  new  satellite  operates  on  what  is  called  an  AA  mode  of  trans¬ 
mission  while  the  existing  McIDAS  equipment  can  only  handle  the  A  mode. 

A  plan  was  developed  to  purchase  a  VAS  frame  synchronizer  and  associated 
circuitry  for  use  with  the  current  data  acquisition  components  in  order  to 
receive  the  new  mode  of  transmission.  Since  the  operational  VAS  satellite 
will  ultimately  replace  the  present  satellite  as  the  primary  data  down¬ 
link,  appropriate  system  modification  is  mandatory. 

^See  Gerlach ,  ed.  (1980),  pp.  SO  -  32,  for  description  of  earlier  acqui¬ 
sitions  . 


Another  modification  was  the  recent  installation  m  two  24-houj 
timers  in  the  Sony  Archive  System.  The  timers  are  attached  via  a  relay 
matrix  to  the  A.C.  [tower  inputs  of  the  two  Sony  2800  V i deo  Recorders. 

The  recorders  can  be  scheduled  to  operate  at  selected  time  periods  by  the 
insertion  or  withdrawal  of  timer  pins.  The  timers  also  have  day-select 
pins  enabliaq  them  to  be  set  up  for  operation  several  days  in  advance. 

This  is  especia.ly  useful  at  night  or  on  weekends  when  an  operator  is  not 
available  yet  satellite  data  acquisition  is  necessary. 

The  Sony  Archive  System  playback  capabilities  were  enhanced  by  the 
addition  of  a  Search  Track  Monitor.  The  monitor  shows  the  position  of  the 
tape  durinq  playback  by  displaying  the  time  of  day,  scan  count,  Julian 
date,  and  year.  This  information  is  recorded  on  a  separate  track  of  the 
video  tape  during  the  real-time  recording  of  the  satellite  transmission. 
When  played  back,  the  track  information  is  sent  to  the  monitor  where  it  is 
decoded  and  displayed  through  the  LED  panel.  The  monitor  is  particularly 
helpful  during  data  extraction  when  the  user  is  looking  for  a  specific 
portion  of  a  tape  which  can  contain  up  to  five  hours  of  recorded  video. 

Documentation  of  McIDAS  hardware  was  begun  in  1980.  Since  its  in¬ 
stallation  in  1974  the  system  has  undergone  numerous  modifications,  includ¬ 
ing  wiring  and  equipment  changes  which  are  not  reflected  in  the  existing 
documentation.  The  intent  is  to  consolidate  in  book  form  documentation  of 
all  the  electronics  in  the  system,  thus  reducing  the  search  time  spent  by 
a  technician  in  isolating  a  malfunction.  Documentation  nearing  completion 
at  contract  end  included  the  data  acquisition  system,  layout  of  the  inter¬ 
connecting  cables  of  the  video  subsystem,  and  construction  of  schematics 
ov  lining  various  control  words  in  the  digital  interface  drawer. 

II.  McIDAS  Support  to  Research  Projects 

In  1980  McIDAS  supported  three  AFGI,  research  projects  using  aircraft: 
the  Large  Scale  Cloud  System  (L.SCS)  study  conducted  by  the  Cloud  Physics 
Branch;  the  Icing  program  conducted  by  the  Climatology  and  Dynamics  Branch; 
and  the  Tactical  Weapons  Support  System  program  conducted  by  the  Optical 
Physics  Division.^ 

2 

See  Gerlach,  ed .  (1980),  pp.  33  -  41,  for  details  of  earlier  support. 
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The  LSCS  program,  a  continuing  AFGL  effort,  sampled  largo  homoqenoous 
cloud  masses  as  they  proqressed  from  west  to  east  across  the  US.  The  air¬ 
craft,  a  specially  instrumented  C-130,  collected  hydrometeor  data  such  as 
particle  size,  number,  concentration,  type,  spatial  distribution,  and  li¬ 
quid  water  content  at  various  altitudes.  Me  IDA:',  satellite  imaqery  loops 
(both  visible  and  infrared)  which  depict  cloud  qrowth  and  decay  as  well  as 
cloud  top  temperatures,  enabled  the  qround  support  team  to  ascertain  the 
most  desirable  locations  for  aircraft  cloud  samplinq.  This  information 
was  relayed  to  the  in-fliqht  team,  allowinq  effective  utilization  of 
available  aircraft  time.  McIDAS  provided  conventional  surface  observa¬ 
tion  plots  of  wind  vectors,  precipitation  amounts,  current  weather,  cloud 
cover,  and  visibility,  along  with  analyses  of  pressure  fields,  streamlines, 
and  dew  point  temperatures.  Upper  air  height  analyses,  radiosonde  plots 
(Stuve' ,  and  cross-section  analyses  also  helped  determine  appropriate 
sampling  altitudes.  Pre-flight  and  in-fliqht  surface  and  upper  air  data 
were  saved  on  magnetic  tape.  Cassettes  containing  satellite  imagery  from 
the  Sony  archive  system  were  saved  for  later  examination.  The  archived 
data  are  then  used  to  document  the  synoptic  situation. 

The  mode  of  operation  was  similar  for  the  Icing  program,  a  study  di¬ 
rected  at  observing  and  measuring  the  accumulation  of  ice  on  an  in-fliqht. 
aircraft.  Using  the  same  aircraft  as  in  the  l.SCS  study,  .'lights  weie  ..-on- 
ducted  to  sample  over  any  one  of  12  radiosonde  stations  within  a  300-mile 
radius  of  Wr ight-Patterson  AFB  (W-PAFB) ,  Oil.  Sequences  of  McIDAS  imagery 
(visible  and  1R)  ,  f>-hour  surface  station  listings,  SLuvo  diagrams,  upper 
air  height  analyses,  and  the  new  icing  program  plots  helped  determine  if 
suitable  icing  conditions  were  present  in  the  study  area.  Again,  pre¬ 
flight  and  in-flight  data  were  archived  on  tape  and  printouts  for  later 
inspection . 

The  Tactical  Weapons  Support  System  program  is  aimed  at.  developing  a 
Tactical  Decision  Aid  (TDA) ,  which  will  use  weather  and  other  geophysical 
information  to  determine  and  predict  if  certain  environmental  factors; 

(e.g.,  1R  target-to-background  contrast,  atmosphei ic  transmission)  jre 
conducive  to  the  deployment  of  specific  precision-guided  munitions  in  the 
tactical  situation.  Using  an  instrumented  KC-115  aircraft  operating  out 
of  W-PAPB,  a  test  flight  was  conducted  in  September  to  obtain  scene  con¬ 
trast  data  over  a  selected  target  area.  McIDAS  was  able  to  provide  visible 
and  infrared  satellite  imagery  loop;;,  .rtaee  data  listings  for  30  stations 
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in  and  adjacent  to  the  target  area,  upper  air  height  analyses,  and  Stuve 
plot  displays  from  the  Huntington,  WV,  Dayton,  OH  and  Pittsburgh  radio¬ 
sonde  stations.  Data  were  collected  for  a  week  prior  to  the  flight  to 
establish  weather  trends.  Appropriate  pre-flight  and  in-flight  data  were 
archived  on  tapes  and  printouts. 

Daily  archiving. of  GOES/SMS  satellite  transmissions  began  in  1980 
with  acquisition  of  the  Sony  cassette  recorders.  By  means  of  two  24 -hour 
timers,  one  recorder  was  set  to  operate  from  1200Z  to  1700Z  and  the  other 
recorder  from  1700Z  to  2200Z.  Each  cassette  tape  holds  five  hours  of  satel¬ 
lite  transmissions.  The  tapes  are  achived  for  approximately  a  month  and 
then  recorded-over  if  there  were  no  requests  for  particular  data.  Tapes 
containing  interesting  weather  phenomena  (hurricanes,  squall  lines,  se¬ 
vere  storms,  volcano  effects)  are  permanently  archived.  Additional  data 
archiving  is  accomplished  with  the  Off-line  Data  Ingest  System  (ODIS) ,  a 
microprocessor-controlled  tape  drive  configured  as  a  peripheral  device  of 
the  McIDAS  facility.  The  ODIS  collects  one-mile  visible,  four-mile  IR 
images  centered  over  eastern  Pennsylvania,  at  hourly  intervals  from  1400Z 
to  2100Z.  Data  are  collected  at  half-hourly  intervals  if  particularly 
interesting  or  severe  weather  is  present  in  the  study  area.  Data  are  then 
analyzed  on  the  AFGL  central  computer  to  calculate  brightness  counts, 
cloud  motion  vectors,  and  IR  temperatures. 

Il  l  .  Syr, tern  Operability 

During  1980  the  McIDAS  facility  was  operated  243  out  of  a  possible 
248  scheduled  working  days,  for  an  availability  of  98  percent  (see  Figure 
D-l) .  Of  the  five  "down"  days,  only  one  was  due  to  a  component  failure; 
the  other  four  were  caused  by  painting  of  the  24-ft.  antenna  dish,  which 
prevented  reception  of  satellite  transmissions. 

IV.  Applications  Software 

SASC  McIDAS  software  efforts  consisted  of  three  types  of  activities: 
debugging  existing  software,  expanding  the  capabilities  of  operational 
programs  t  lirougli  add  it  ion  of  new  options,  and  development  of  new  proqrams. 

The  main  body  of  system  and  applications  software  for  the  McIDAS  was 
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written  at  the  University  of  Wisconsin,  whose  Space  Science  and  Engineer¬ 
ing  Center  developed  the  McIDAS.  However,  many  modules  required  modifi¬ 
cation  to  allow  for  configurational  differences  at  AFGL.  Others  contained 
errors  introduced  by  imperfect  copying  of  the  original  proqrams,  and  a  few 
of  the  most  recent  additions,  received  from  the  University  of  Wisconsin  in 
an  early  stage  of  development,  required  correction  of  logic  faults. 

In  order  to  adapt  the  original  applications  software  to  the  special¬ 
ized  needs  of  AFGL  and  SASC  scientists,  new  options  were  added  by  SASC 
to  many  existing  programs.  These  additions  increased  the  versatility  of 
NOWCASTING  proqrams  (used  to  decode,  store,  analyze,  and  display  conven¬ 
tional  surface  and  upper  air  data  received  via  a  high  speed  FAA  data  line) , 
expanded  the  graphics  capabilities  to  include  storage  on  television  frames 
through  a  plot  display  file,  and  improved  speed  and  accuracy  of  the  navi¬ 
gation  fine-tuning  process. 

New  programs  unique  to  the  AFGL  McIDAS  and  specifically  designed  to 
support  current  research  projects  were  incorporated  into  the  system.  In 
order  to  expand  the  data  base  upon  which  imagery  information  extraction 
techniques  are  tested,  data  ingestion  software  for  two  geostationary  satel¬ 
lite  imaqery  archives  (a  30-day  rotating  archive  of  full  resolution  hemi¬ 
spheric  data  from  the  video  cassette  archiving  system,  and  a  permanent 
ongoing  collection  of  9-track  tapes  covering  the  New  England  region  re¬ 
corded  by  the  ODIS)  was  developed. 

A  new  program  by  SASC  was  added  to  McIDAS  to  assist  forecasters  in 
estimating  the  probability,  type,  and  intensity  of  aircraft  icing  at  vary¬ 
ing  altitudes.  The  algorithm  developed  was  based  on  a  sequence  of  opera¬ 
tional  procedures  used  by  Air  Weather  Service  forecasters . ^  Upper  air 
sounding  data  ate  used  to  estimate  liquid  water  content,  layer  stability, 
and  temperature  for  correlation  with  empirically  derived  icing  probabili¬ 
ties  . 

A  major  addition  to  McIDAS  software  was  the  Mesoscale  Forecasting 
Facility  (MFF) ,  developed  by  SASC  as  part  of  a  program  of  preliminary 
testing  and  examination  of  the  feasibility  of  using  a  man-interactive  sys¬ 
tem  in  the  preparation  of  short-range  terminal  forecasts.  McIDAS  was  used 
to  simulate  many  aspects  of  the  interactive  mi nicomputer-based  display  sys- 

AWS ,  1980:  Forecasters '  guide  on  aircraft  icing.  Technical  Report 
AWS/TR-80/001 .  Scott  AFB,  II,. 
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tem  envisioned  for  the  Automated  Weather  Distribution  System  (AWDS) 
planned  for  Air  Weather  Service's  Base  Weather  Stations  by  the  mid- 
1980' s  and  the  Satellite  Data  Handling  System  (SDHS)  planned  for  the 
Air  Force  Global  Weather  Central. 

The  MFF  is  a  "menu-driven"  system  designed  to  assist  in  preparing, 
executing,  and  verifying  short-range  forecasts.  A  series  of  questions  and 
instruction  pages  is  generated  which  assists  the  forecaster  in  examination  of 
current  weather  data,  formulation  of  1-,  3-,  and  6-hour  probability  fore¬ 
casts,  insertion  of  predictions  into  personalized  data  files,  verification 

4 

of  forecasts,  and  later  recall  for  evaluation  ol  performance. 

A  nine-month  test  ( January-September  1980)  was  conducted  by  AFGI,  and 
SASC  meteorologists  to  evaluate  the  utility  of  the  interactive  approach. 
Fifteen  researchers  with  varying  degrees  of  forecasting  expertise  partici¬ 
pated  in  the  program.  Three  teams  performed  the  forecasting  functions  on 
a  weekly  rotating  basis. 

The  exercise  was  structured  in  the  form  of  a  daily  map  discussion 
led  by  the  Chief  Forecaster  of  the  Day,  a  role  assumed  by  each  team  member 
one  day  during  the  week.  The  Chief's  primary  resionsibil ity  was  prepara¬ 
tion  of  a  forecast  problem  for  consideration  by  the  team.  Based  on  an 
assessment  of  that  day's  weather,  the  Chief  selected  a  station  which  pro¬ 
vided  an  interesting  mesoscale  forecasting  situation,  and  three  forecast 
components.  Predictand  options  included  wind  speed,  wind  direction,  gusts, 
visibility,  cloud  cover,  ceiling  height,  temperature,  dew  point  tempera¬ 
ture,  and  weather.  A  map  file  of  appropriate  weather  information  was  pre¬ 
pared  for  use  in  the  group  discussion.  Data  sources  available  through  the 
MFF  included  conventional  surface  and  upper  air  observations  (Services  A 
and  C)  and  real-time  GOES  satellite  imagery  (visible  and  infrared  channels). 

A  wide  range  of  objective  analysis  and  display  routines  was  provided  alonq 
with  image  looping  and  data  overlay  capabilities  for  combining  convention¬ 
al  and  satellite  data.  Following  the  map  discussion,  each  forecaster 
entered  predicted  values  of  the  assigned  parameters  for  forecast  lengths 
of  1,  3,  and  6  hours.  Verification  and  scoring  of  individual  forecasts 
were  performed  the  following  day. 

At  the  end  of  the  test  period,  a  constructive  review  of  the  experi¬ 
ment  was  conducted  to  evaluate  the  advantages  and  disadvantages  of  an  AWDS/ 

4 

See  Gerlach,  ed.  (1980),  pp.  49-73,  f  a  detailed  description  of  the  MFF. 
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SPHS  Ml-T  system  !  'i  i  >i  e  i  .1 1  i ,  .11.1 1  met -Wiit  eli  and  short-range  prediction  .  An 
.issi'ssmt'iu  w.is  also  modi'  . .!  intuit*  direct  ions  the  protect  should  take 
tow. ltd  the  oli  irrt  i  ves  oi  t  ■  1 1  mu  l  .it  i  nq  and  testing  mesoscale  nowcasting  and 
toioiMst  i  ns  techniques  tor  the  AWDS/SDUS  era. 

Among  'he  advantages  cited  by  MFI-*  participants  was  the  flexibility  is. 
analyzing  aid  display  in>i  conventional  data  made  possible  through  an  inter¬ 
act  ive  s.yst  -in.  ilouily  surt  a--e  cl. aits.,  regional  dis(  lays,  act  .a'  e.i  •  so, ms- 
ings,  and  d  a  ivea  tields  s.ucn  as  divergence,  vorticity,  and  temperature 
adve  t ion  are  examples  of  the  tailored  products  shown  to  be  useful  that 
are  sot  available  on  standard  facsimile  charts.  Tnstant  recall  of  observa¬ 
tion'.  ordered  by  time,  station,  region,  or  weather  type  was  found  to  be  an 
improvement  over  conventional  modes  of  display.  The  addition  of  real-time 
satellite  imagery  alone,  in  time  lapse  sequences,  and  in  combination  with 
surf  ice  and  RAOB  reports,  was  a  valuable  expansion  of  the  mesoscale  data 
base.  Other  comments  noted  the  advantages  of  the  plot  display  file  for 
instantaneous  recall  of  displays  and  time  sequencing  of  surface  observa¬ 
tion  i,  and  reduced  or  eliminated  dependence  on  paper  products. 

Disadvantages  included  the  deficiency  in  mesoscale  information 
caused  by  coarse  spacing  of  the  operational  reporting  network  and  lack  of 
radar  data.  Bad  navigation  on  rare  occasions  rendered  satellite  imagery 
misleading.  Several  of  the  more  experienced  forecasters  had  problems  ad¬ 
just  i.ng  to  the  new  display  format,  expressing  difficulty  in  obtaining  an 
inti  (rated  understanding  of  the  weather  situation  when  limited  to  one  chart 
at  ci  time.  It  was  suggested  that  the  addition  of  one  or  more  monitors 
woull  alleviate  the  problem  by  allowing  coincident  display  of  several  analy¬ 
ses.  Lin  rations  in  screen  size,  resolution,  and  color  options  were  men- 
tior  id  ai  undesirable  features  of  the  system,  along  with  such  mechanical 
difficulties  as  occasional  flooding  of  the  screen  and  delays  while  graphics 
were  produced.  It  was  suggested  that  a  means  of  subjectively  modifying 
objc  rtively  analyzed  fields  was  needed  to  take  advantage  of  independent  in- 
fornation  and  the  experience  of  the  forecaster. 

Examination  of  verification  statistics  for  the  test  period  revealed 
a  g<  leral  trend  toward  decreasing  skill  with  increasing  forecast  length, 
as  expected.  A  wide  range  of  predictability  amonq  the  various  forecast 
parameters  was  indicated,  with  greatest  success  achieved  in  predicting  dew 
!x>int  temperature,  visibility,  and  temperature,  and  least  skill  in  predict¬ 
ing  ceiling  height  and  wind  direction.  Firm  conclusions  were  restricted 
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by  the  small  sample  size.  j 

The  MFF  experience  was  a  successful  first  step  in  evaluating  the 
advantages  and  disadvantages  of  an  interactive  facility  for  producing 
short-range  forecasts.  The  resulting  recommendations  provide  a  basis  for 
continued  research  into  the  kinds  of  data,  displays,  and  capabilities  which 
will  be  supplied  to  AWDS-era  forecasters. 


i 

I 

I 

i 


51 


E.  WEATHER  RADAR  TECHNIQUES 


I.  Wind  Retrieval 

Ability  to  lerivo  automatically  the  fields  of  motion  within  and 
about  severe  stc  rms  during  sinqle  radar  operations  is  highly  desirable. 
Several  approacl es  to  this  problem  have  been  proposed,  but  as  yet  none 
lias  been  tested  operationally. 

Smythe  (19110)^  presented  a  pattern  recognition  technique  for  re¬ 
trieving  nearly  horizontal  velocity  vectors  in  clear  air  by  cross- 
con  alating  data  from  two  consecutive  scans.  Even  though  the  Doppler 
radar  is  blind  to  motions  perpendicular  to  its  narrow  beam,  when  this 
technique  is  used,  cross-radial  as  well  as  along-radial  components  of 
velc  :ity  may  be  derived. 

Software  associated  with  the  wind  retrieval  technique  was  modified 
for  and  coded  on  the  Interdata  7/12  minicomputer  system  at  the  AFGL 
Weather  Radar  Facility.  In  the  process  of  implementing  this  software, 
several  problem.*  were  encountered  by  SASC  radar  meteorologists.  These 
were  in  part  dui  to  system  size  and  compiler  differences. 

AJ.  10 ,  durii  q  implementation  the  software  was  modified  to  improve 
efficiency.  FailLer,  three  jobs  were  run  to  obtain  the  finished  product; 
at  AFGL,  the  first  two  programs  were  merged  by  SASC  and  the  process  was 
reduced  to  two  jobs.  One  proqram  prepares  digitized  data  for  the  second 
program's  correlation  analysis  and  transrorms  the  digitized  data  to 
scalar  data. 

Two  software  versions  of  the  technique  had  been  previously  developed. 
One  correlates  DoppLer  t.iiia]  velocities,  the  other  reflectivities  or 
spectral  widths.  !n  both  versions,  radial  velocities  are  used  to  deter¬ 
mine  radial  displacements  of  small  regions  during  the  time  elapsed  between 
scans.  Modification  ot  the  first  proqram,  which  transforms  the  digitized 
da'  a  to  scalar  data,  was  initiated  by  SASC  to  provide  the  analysis  proqram 
with  either  velocities  and  reflectivities  or  velocities  and  spectral 
widths.  When  finished,  the  modified  proqram  could  bo  interfaced  with  the 

'smythe,  G .  R .  ,  l'JHi):  The  retrieval  of  mososcale  horizontal  winds  using  a 
sinqle  Doppler  radar.  I’reprjnts,  19th  Conference  on  Radar  Meteorology, 
AMS,  Boston;  pp.  250-255. 
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reflectivity  and  spectral  width  correlation  program.  Wind  fie  Is  c<uld 
then  be  produced  from  the  ensuing  analyses. 

The  pattern  recognition  technique  was  applied  to  data  col  ecte<  on 
May  2,  1978  in  Norman,  OK  during  the  approach  and  passage  of  a  gust  front. 
Results  from  several  analyses  indicated  the  presence  of  the  gu  ;t  fr<  nt . 
Doppler  radial  velocities  were  correlated  for  these  analyses. 

Development  of  a  real-time  version  of  the  technique,  regu  ring  less 
core  memory,  was  also  initiated  by  SASC.  This  involves  dividi  ig  tht  analy¬ 
sis  program  into  two  parts:  the  first  program  partitions  the  first  scan  in¬ 
to  small  regions,  computes  the  variance  and  radial  displacemen:  of  each  re¬ 
gion,  and  stores  this  information  for  each  small  region  on  sep  irate  disk 
files.  This  processing  is  completed  before  the  second  scan  is  colli cted. 

The  second  program  uses  the  stored  information  to  find  equally  di¬ 
mensioned  regions  of  the  second  scan  (collected  at  the  same  elevation  angle 
as  the  first  scan)  and  determines  the  covariances  between  data  of  t)  ese 
regions  and  data  of  radially  and  azimuthally  displaced  areas  of  the  first 
scan.  It  also  computes  variances  of  the  second  scan  regions  a  ad  soi ts  and 
evaluates  the  many  correlation  coefficients  computed.  Finally,  it  pro¬ 
duces  a  display  of  the  vector  wind  field. 

Additionally,  at  the  end  of  the  contract  period  two  curves  were 
being  least-squares  fitted  to  the  correlation  data.  The  curves  are  an 
exponentially  decaying  cosinusoid  and  a  third  order  polynomial.  It  was 
expected  that  one  curve  or  the  other  would  properly  smooth  the  correlation 
field  and  perhaps  give  a  more  precise  estimate  of  the  locations  of  t he  end 
points  of  the  wind  vectors. 

II.  Real-Time  Storm  Cell  Detection  and  Tracking 

This  SASC  effort  was  to  develop  a  real-time  software  technique  for 

the  detection  and  tracking  of  significant  storm  cells.  The  processing 

2 

algorithm  originally  developed  by  Crane"  extracts  the  essential  information 
required  to  describe  the  organization,  development ,  and  motion  of  tie  storm 
system.  This  information  includes  the  characteristics  of  the  uctiv.-  cells 

2 

Crane,  R.K.,  1979:  Automatic  cell  detection  and  trackmq,  IEb'E  Trvns. 
on  Geoscience  electronics,  (IE- 17,  250-2C2. 
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within  tiie  storm  system,  the  clusters  of  cells,  the  fixed-threshold  con¬ 
tours  or  echo  regions  that  encompass  the  cells,  and  the  complexes  of  echo 
regions  that  comprise  the  storm  system. 

The  algorithm  was  modified  in  1079  for  use  on  the  Interoata  7/37  on¬ 
site  com  iuter  of  AFGL.  This  modified  (1979)  program  is  aide  to  opei.ro 
over  a  u ide  ranqe  of  reflectivity  thresholds.  In  order,  however,  to  re¬ 
duce  tht  number  of  cells  processed  in  the  computer,  thereby  imj  r<-v;  r.  : 
operatir  speed  for  real-time  operations,  SASC  raised  the  reflect! vi: v 
threshol  to  40  dRz.  Despite  this  improvement,  the  current  running  time 
for  refl  ctivity  processing  is  130  seconds  per  azimuth  scan.  This  is  not 
satisfac  ory  for  real-time  operation,  which  requires  40  seconds  per  azi¬ 
muth  sea;  . 

The  SASC  developmental  work  was  therefore  organized  into  three  tasks: 

a.  To  check  the  accuracy  of  the  original  algorithm, 

b.  To  implement  the  modified  (1979)  program  for  real-time 

use,  and 

c.  To  test  yet-to-be-developed  software  for  real-time  opera¬ 
tion  . 

Under  tb  contract  efforts  were  directed  toward  tasks  a  and  b. 

The  first  step  in  task  a  was  to  provide  data  from  the  pulse  pair  re¬ 
corder  i  terface  (PPRI)  as  input  to  the  original  tracking  algorithm,  us- 
inq  a  pr  gram  that  transforms  digitized  data  to  reflectivity  values.3 
This  dig  tized  program  is  also  capable  of  decoding  the  variance  and  velo¬ 
city  va!  es .  A  double  buffering  technique  developed  by  SASC  was  intro¬ 
duced  in  o  the  digitized  program  to  increase  input/output  speed.  A  mag¬ 
netic  di  ;k  is  used  as  the  output  device  because  the  data  transfer  speed  of 
the  tape  drive  ;  relatively  slow.  Therefore,  two  additional  programs  were 
developed:  DSCTOTP,  which  transfers  the  data  from  disk  to  tape,  and  DIGRD, 
which  prints  the  data  from  disk  or  tape. 

Concurrently,  the  output  format  of  the  digitized  program  was  modified 
to  conform  to  the  Common  Doppler  Radar  Exchange  Format  adopted  by  the 
Doppler  radar  community.  This  involved  adding  seven  ancillary  data  header 
blocks  to  the  decoded  radar  data  blocks. 

For  task  b,  the  logic  of  the  modified  (1979)  program  was  examined  and 

3Roak,  T.I.S.,  /...I.  Jagodnik,  R.B.  Marshall,  D.  Riceman,  and  M.  Young, 

1977:  Tracking  .md  significance  estimator.  Final  Ret>ort,  Raytheon  Company. 
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redundant  execution  steps  were  removed.  Timinq  analysis  showed  that  the 
speed  of  the  data  acquisition  subroutine  and  peak  calculation  subroutine 
needed  to  be  increased  for  real-time  implementation.  Since  codinq  these 
subroutines  in  assembler  language  would  reduce  the  time  required  for  in¬ 
direct  storage  addressing  and  bit  manipulation,  SASC  developed  Common 
Assembly  Language  (CAL)  macros  for  them.  Such  program  modularization, 
that  is,  dividing  large  programs  into  subprograms  which  are  developed  and 
debugged  individually,  conforms  to  good  programming  pract  ice. 

III.  Detection  of  Turbulence  in  Thunderstorms  using  Radar  Methods 

(A) .  Introduction 

The  detection  of  turbulence  in  thunderstorms  by  remote  methods  has 
been  a  subject  of  investigation  for  many  years.  Early  studies,  relying 
on  observations  from  penetrating  aircraft,  attempted  to  correlate  turbu¬ 
lence  encounters  with  radar  storm  features.  The  general  conclusion  was 
that  the  frequency  of  occurrence  of  moderate  to  severe  turbulence  increased 
markedly  with  increasing  storm  radar  echo  intensity  (dBz) ,  with  40  dBz  a 
rough  threshold  value  for  development  of  such  turbulence  somewhere  in  the 
storm.  Further  observation  indicated  that  severe  turbulence  was  widely 
distributed  in  storms  whenever  the  intensity  indicated  damaging  hail 
(about  50  dBz) ,  with  turbulence  encounters  outside  high  intensity  cores 
being  as  frequent  as  those  within  (Lee,  1967).^  This  was  particularly 
evident  within  thunderstorm  lines  where  severe  turbulence  was  frequently 
encountered  in  low  echo  intensity  regions  between  individual  storms.  The 
conclusion  was  that  storm  echo  intensity  was  an  indicator  of  the  presence 
of  turbulence,  but  it  did  not  reliably  measure  its  strength  nor  locate  its 
position. 

In  light  of  these  results  and  with  continued  development  in  Doppler 

radar  techniques,  emphasis  was  placed  on  finding  kinematic  signatures  of 

2 

hazardous  zones  within  thunderstorms,  Armstrong  and  Donaldson  (1969) 
developed  the  Plan  Shear  Indicator  (i'SI),  which  was  able  to  detect  regions 

^Lee,  J.T.,  1967:  Atmospheric  turbulence  and  radar  echoes,  in  Technical 
Memorandum  IERTM  -  NSSL,  32^  Norman,  OK,  pp.  1-9. 
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Armstrong,  G. ,  and  R.  Donaldson,  dr.,  1969:  Plan  Shear  Indicator  for  real¬ 
time  Doppler  radar  identification  of  h  .zardous  storm  winds,  d .  Appl .  Meteor . , 
8,  pp.  176-18  1. 


<>t  st  roiv,  MhtMr  of  I  ht'  India]  wind.  Leo  .uid  Kraus  (1975)  \  by  comparing 
I'SI  and  penetrating  aircraft  observations,  found  good  correlation  between 
reqions  of  strong  sliear  of  the  radial  wind  and  severe  turbulent  events. 

Radar  data  of  Frisch  and  Strauch  (1970)  suggested  that  hiqh  turbulent  inten¬ 
sity  may  not  necessarily  be  correlated  with  zones  of  hiqh  shear  of  the  ra¬ 
dial  wind,  but  when  present,  may  be  found  in  reqions  exhibitinq  larqe 

4 

Doppler  spectrum  variance.  Later  Lee  (1977)  also  found  qood  correlation 
between  reqions  of  severe  turbulence  and  larqe  spectrum  variance.  Thus, 
present  attempts  arc  beinq  directed  at  usinq  Doppler  spectrum  variance  as 
a  turbulence  indicator. 

The  turbulent  contribution  to  Doppler  spectrum  variance  is  a  function 
of  the  turbulent  intensity  (or  alternatively  the  eddy  dissipation  rate  and 
the  turbulent  outer  scale) ,  range  of  the  observation  volume  from  the  radar, 
and  precipitation  environment.  The  manner  in  which  these  effects  combine 
is  not  well  known,  and  thus,  the  interpretation  of  this  indicator  is  not 
yet  completely  understood.  The  aim  of  this  research  is  to  clarify  the 
relationship  betwee  i  Doppler  spectrum  variance  and  these  parameters  and  to 
develop  a  more  reliable  radar  turbulence  detection  technique. 

This  report  focuses  on  results  of  theoretical  investiqations  and  data 
analyses  performed  oy  SASC  durinq  the  period  May,  1979  to  June,  1980. 

Primary  topics  discussed  are:  l)  the  response  of  a  distribution  of  precipi¬ 
tation  particles  to  various  ranges  of  scales  of  turbulent  motion;  2)  the 
relationship  between  eddy  dissipation  rate  and  Doppler  spectrum  variance 
as  a  function  of  the  turbulent  scale  regime  and  precipitation  environment; 
and  3)  the  distribution  of  local  turbulence  regimes  within  thunderstorms 
as  Ue Lineated  by  aircraft  data.  Finally,  since  one  of  the  ultimate  goals 
of  this  work  is  to  levise  a  methodology  for  turbulence  detection  by  inco¬ 
herent  radar,  a  [xnssible  method  by  which  non-Doppler  radar  may  estimate 
wind  shear  for  the  purpose  of  removing  its  contaminating  effects  from  the 
Doppler  spectrum  variance  will  be  briefly  discussed.  Results  of  these  in- 

’  I.re,  J.T.,  and  M.  Kraus,  1975:  Plan  Shear  Indicator  and  aircraft  measure¬ 
rs  its  of  thunderstorm  turbulence:  experimental  results.  Preprints ,  16th 
Radar  Meteorology  Conference,  AMS,  Boston;  pp.  337-340. 
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Lee,  J.T.,  1977:  Application  of  Doppler  weather  radar  to  turbulence 
measurements  which  affect  aircraft.  Final  Report  FAA-RD-77-145  to  Systems 
Research  and  Development  Service,  FAA ,  Washington,  D.C.,  20  pp. 
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vestigations  point  to  the  need  for  further  data  analysis  for  verification 
and  selection  of  practical  solutions. 


(B) .  Theoretical  Investigations 


The  response  of  precipitation  particles  to  turbulent  air  motions  is 
of  fundamental  importance  since  the  Doppler  radar  velocity  is  generally 
equated  to  the  radial  component  (component  along  the  radar  viewing  direc¬ 
tion)  of  the  environmental  wind.  However,  the  Doppler  velocity  is  actually 
the  power-weighted  mean  of  the  Doppler  spectrum,  which  is  formed  from  mea¬ 
surements  of  the  radial  velocities  of  the  precipitation  particles  within 
the  radar  pulse  volume,  with  the  spectral  power  at  each  spectral  velocity 
value  proportional  to  the  sum  of  the  sixth  powers  of  the  diameters  of  the 
precipitation  tracers  moving  with  that  radial  velocity  value.  It  is  sus¬ 
pected  that  a  sizeable  portion  of  Doppler  spectrum  power  may  be  supplied  by 
tracers  which,  because  of  their  large  inertia  (e.g.,  large  raindrops),  do 
not  follow  turbulent  fluctuations  faithfully.  If  so,  then  an  incorrect 
estimate  of  the  turbulent  motion  (e.g.,  Doppler  spectrum  mean  and  variance) 
is  obtained.  Thus,  one  is  interested  in  determining  the  resjxanse  of  a  dis¬ 
tribution  of  particles  to  various  scales  of  motion.  This  may  be  obtained 
by  appropriately  summing  the  individual  responses  of  the  particles  com¬ 
prising  the  distribution. 

Consider  a  single  particle.  Its  equation  of  motion  may  be  written  as 


dV  dV  p  +  C  N 

m  — =  m  +  m  g(l  -  —  )  -  Gnap(V  -  V  )  — 

D  dT  a  dT  D  D  a  2 4 


(1) 


where  m.,  V,,  p.,  and  a  are  the  mass,  velocity,  density,  and  radius  of  the 
d  d  d  ^ 

particle.  Air  density  is  p  ,  m  is  the  mass  of  displaced  air,  and  V  is 

a  a  a 

the  air  velocity  in  the  vicinity  of  the  particle.  Kinematic  viscosity  is 

— ► 

U,  gravitational  acceleration  is  g,  and  and  (=  2  Re)  are  the  drag 

coefficient  and  alternate  Reynolds  number  of  the  particle.  The  coefficient 

C,N  /24  is  a  catch-all  term  and  is  actually  related  to  fluid  uccelerut ion 
d  Re 

and  past  particle  motion  (Pruppacher  and  Klett,  1978).  This  equation 
simply  relates  the  net  force  acting  on  the  particle  to  the  pressure, 
buoyancy-corrected  gravity,  and  drag  terms.  Equation  1  may  be  written  as 
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Iruppaehor ,  II. R.,  and  ,1.1).  Klett,  197  Mi  rophysic:.  of  Clouds  and  Pre¬ 
cipitation,  Huston,  I).  Keidel  I  'ill  >  I  i :  :|i  i  I "  I  Co.,  VIS  pp. 
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Acti  ,j1  solution  of  liquation  2  requires  appl  Lea t ion  of  numerical 

>  > 

methods,  since  C  is  dependent  on  (V.,  -  V  ).  To  obtain  <\  qencr  •. i  -;u  jyr  : 

i  J  l  >  d 

solution  for  the  particle  rc-s;onsr  it  i  necessary  to  make  some  sin.;  ii  ty¬ 
ing  assumption  concerninq  the  draq  force.  For  this  purpose  it  is  assumed 
that,  on  the  averaqe,  the  particle  is  settling  at  its  terminal  fallspeed 
(V  ) *  The  drag  force  in  Equation  2  may  be  replaced  by  its  average  value. 
Then  following  Stackpole  (1961)  the  following  relation  is  applicable: 


Neglecting  particle  growth  and  the  variation  of  air  density  with  height, 

T  is  a  constant  for  each  particle  .  Note  that  this  assumption  means  that 
the  draq  coefficient  is  replaced  by  the  value  attained  when  the  particle  is 
settling  at  its  normal  terminal  fallspeed  in  still  air,  and  thus  it  is  a 
slightly  averaged  particle  response  that  may  be  estimated. 

In  .ipplyinq  Equation  2  to  a  turbulent  field,  one  generally  considers 
turbulence  as  composed  of  a  random  collection  of  eddies  having  a  range  of 
scales  (o.g.,  wavelengths  X).  Mathematically,  one  may  view  this  collec¬ 
tion  as  a  series  of  harmonic  functions  of  varying  frequency  (to)  or  wave- 
number  (K  =  2ir/A) ,  having  random  amplitude  and  combining  to  produce  the 
observed  turbulent  structure.  The  response  of  a  particle  to  a  turbulent 
gust  may  therefore  be  considered  the  result  of  a  superposition  of  the 
Pa.  tide's  responses  to  this  collection  of  harmonic  functions,  or  eddies. 

The  form  of  the  particle  response  to  a  single  such  function  may  be 
seen  by  considering  one  dimension  (horizontal)  and  introducing  an  air 

velocity  V  =  V  cos  (OT  (5) 

a  o 

into  Equation  2.  The  f .article  velocity  is  found  to  be 

T 

2  2  ?  ^  2  2  ^  _ 

Vd  =  V  [(]  t  <0  T  f/)  V  (l  +  (0  T  ) '  J  cos  (WT  +  U)T(P0-  D)  +  Ae”x  (6) 


fjtackfjole,  J.D.,  1961:  The  effectiveness  of  raindrops 
sors.  Proceedings,  9th  Weather  Radar  Conference,  AMs 
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a 

where  Pq  -  —  and  urr  <  1  is  assumed.  This  last  condition  is  met  for  all 
environments  except  those  combining  lurqe  particle  size  and  turbulent  scale 
<  60  m.  (Even  then,  the  energy  content  of  the  turbulent  field  at  those 
scales  is  small  in  relation  to  larger  scales  and  need  not  be  considered 
further.)  The  first  term  shows  that  the  particle  velocity  is  reduced  in 
amplitude  and  phase-shifted  relative  to  the  air  velocity.  The  second  term 
is  a  transient  term.  (These  effects  are  generally  small.)  In  reality, 
there  is  no  way  of  actually  incor (orating  the  phase  and  transient  informa¬ 
tion  into  actual  motion  of  the  particles  which  enter  any  given  region  of 
turbulence,  since  their  present  motion  is  also  dependent  ujon  their  motion 
prior  to  entering  this  region,  and  this  is  unknown.  Thus,  only  the  ampli¬ 
tude  response  will  be  considered. 

It  is  necessary  to  map  this  response  function  from  dependence  upon 
frequency  to  turbulent  scale,  since  frequency  is  dependent  upon  the  rela¬ 
tive  speed  of  the  atmospheric  structure  and  the  sensing  instrument  (here 
the  particles),  and  therefore  is  an  artificial  turbulent  parameter.  In  a 
manner  similar  to  Taylor's  hypothesis,  one  may  form  a  relationship  between 
turbulent  wavenumber,  radian  frequency,  and  particle  terminal  fallspeed  as 


resulting  in  the  response  function  taking  the  form 


VVa  “  1/(1  +  Kz2vT4/g2)45  (8) 

.  -3 

where  Pq  «  10  has  been  assumed.  This  relationship  was  first  obtained  by 
Stackpole  and  states  that  as  the  particle  size  increases  (V  increases) 
the  response  to  turbulent  motions  having  any  qiven  apparent  scale  A 
(=  K^/2tt)  decreases.  Conversely,  a  particle  of  given  size  (V  constant) 
exhibits  decreasing  response  to  decreasing  scale.  The  term  "apparent 
scale"  has  been  used  to  signify  that  the  particles  are  responding  to  the 
field  which  is  the  projection  of  the  3-D  turbulent  field  upon  the  z-axis 
and  that  there  may  not  be  a  one  to  one  correspondence  of  features  in  the 
3-D  field  to  those  observed  in  the  1-D  field. 

It  is  instructive  to  consider  the  effect  this  particle  response  has 
upon  the  actual  mean  total  velocity  variance  contained  by  a  particle  dis¬ 
tribution  when  in  a  turbulent  field  .  This  is  given  by 


50 


K 

2  2  2  /*  rn^ix 

VAR  =<u  +  v  +w>=/  E(K)W(K)dK  (9) 

Kmin 

where  the  air  enerqy  density  (9)  is  assumed  given  by  the  Kolmoqorov 
relation 

2/3  -5/3 

E(K)  =  Ce  7  K  7  (10) 

in  which  C  is  a  known  constant,  c  is  the  eddy  dissipation  rate,  and 
K  =  | K | .  The  particle  distribution  response  is 


l(K)  =  2J  j  [{qN(D)atan(KV^/q)]/KV^)|  dD 


N (D) dD 


where  the  particle  number  density  and  fallspeed  are  given  by 

N (D)  =  N  e~^D 
o 

and 


(12) 

(13) 


respectively.  Figure  E-l  shows  this  variance  (VAR  ) ,  normalized  to 

pD 


the  air  turbulent  variance  VAR  where 

a 


2 1' 


E (K) dK. 


The  constant  parameters  in  Equations  11-14  are  a  =  1690  cm'4  sec  1 , 

b  -  .6,  b  .  =  .005  cm,  D  =  .5  cm,  and  K  =  2lTcm  ^  (corresponding 

nun  max  max 

to  a  minimum  turbulent  scale  of  1  cm)  .  Parameters  K  .  (=  2TT/X  ) 

min  o 

corresponding  to  a  turbulent  outer  scale  A  and  the  particle  distribution 

o 

slope  A  are  the  independent  variables. 

This  t 1  jure  chows  that  the  mean  turbulent  variance  of  the  particles 
is  typically  95  percent  (or  greater)  of  the  air  value,  for  all  precipita¬ 
tion  environments  (e.g.,  heavy  rain  A  =  10  cm  *  to  light  rain  A  =  60  cm  ^ ) 

and  outer  scales  (A  )  greater  than  500  m. 
o 

This  figure  demonstrates  that  even  though  individual  particles  (e.g., 
large  drops)  may  not  respond  well  to  small  scales  of  motion,  the  actual 
me..  turbulent  velocity  variance  of  the  [or tide  distribution  is  generally 
a  good  measure  of  the  atmospheric  value.  Unfortunately,  the  radar  does 
not  weight  the  particles  identically,  as  this  calculation  assumes,  but 
rattier  weights  tiie  particles  according  to  the  sixth  [»wer  of  drop  diameter. 
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OUTER  SCALE  Lo(Km) 


The  total  3-D  radar  measurable  turbulent  variance  is  given  by 


where 


2  2  2  /*  rnax 

VAR  =  <u^  +  v  +  w^  >  =  /  E(K)W  (K)dK 

P  •'K  .  1 

min 

W2(K)  =  2 ^  j[(qN(D)D6atan(KV^/q)]/KV,^  |  diy/^N 


(15) 


(D) D  dD .  (16) 


This  quantity,  normalized  to  the  actual  air  turbulent  variance,  is 

shown  in  Figure  E-2 .  It  is  observed  that  if  the  outer  scale  (A  )  extends 

o 

to  a  kilometer  or  more,  then  the  radar  can  measure  90  to  97  percent  of 
the  air  motion  variance  in  all  precipitation  environments.  For  shorter 
outer  scales,  however,  particularly  in  heavy  rain  situations,  the  radar 
estimate  may  be  in  significant  error.  As  an  example,  for  A  =  10  cm  ^ 
and  an  outer  scale  of  250  m,  the  total  radar  measurable  variance  is  only 
about  72  percent  of  the  air  value.  The  striking  difference  between  Fig¬ 
ures  E-l  and  F-2  results  from  the  strong  power  contribution  to  the  Doppler 
spectrum  from  large  particles  which  exhibit  poor  response  characteristics, 
and  demonstrates  that  the  usual  assumption  of  precipitation  being  perfect 
tracers  of  the  air  motion  can  lead  to  radar  turbulent  energy  estimates 
having  significant  error.  The  results  demonstrate  that  imperfect  particle 
resfonso  and  th<  finite  scale  of  turbulence  must  be  taken  into  account  if 
reasonable  estimates  of  turbulent  parameters  are  to  be  obtained  by  radar 
methods.  This  will  be  accomplished  by  modifying  the  pertinent  basic  radar 
relations. 

In  all  that  follows  the  wind  environment  is  considered  fully  turbulent 
to  clearly  obtain  the  relationship  between  the  radar  observed  parameters 
and  the  turbulei  t  field.  The  results  are  also  applicable,  however,  in  the 
presence  of  a  constant,  background  wind. 

In  general  form  the  Doppler  spectrum  variance  may  be  written  as 


2  —  2 
VAR  =  V  -  V 
P  P 


(17) 


where  V 

f’ 

over  the 
variance 
square  o 
equation 


is  par'  icle  radial  velocity  and  the  overbar  indicates  averaging 
Dopplei  tower  spectrum.  This  relation  states  that  the  spectrum 
equals  the  averaged  square  of  particle  radial  velocity  less  the 
f  the  u-'eraqed  particle  radial  velocity.  More  specifically,  this 
may  he  written  as 
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OUTER  SCALE  Xo(Km) 


magnitude  (bottom)  and  outer 


VAR  = 


3  JJ V^RJqtV^.R)  I(r,R)dV^dR 
C  j  JJ '  i  ( V  ^  ,  R )  I  (r,R)dV  dR 

C,  //  V  (R)  q  (V  ,  R)  X  (r ,  R)  dV  dR 
!>  P  I> 

cjj  q(V{/R)  i  (?,R)  dv  dR 


showing  more  clearly  the  |>ower-weiqhted  averaqinq  and  pulse  volume 

averaging  over  available  particle  radial  velocities.  Here  C  is  the 

a-  ^ 

radar  constant,  q(V^,R)  is  the  reflectivity  of  the  particles  in  the  vol¬ 
ume  element  dR,  (proportional  to  the  sixth  power  of  the  partial  diameter) , 

-V  ->• 

located  R  from  the  pulse  volume  center  located  at  r  (relative  to  the  radar), 

and  havinq  radial  velocities  ranqinq  from  V  to  V  +  dV  .  The  term 
a  a  P  P  P 

C^q(V  ,R)dRdV ^  is  the  power  back-scattered  by  these  particles  to  the 

radar,  and  I(r,R)  is  the  two-way  beam  illumination  pattern  which  is  as- 

sumed  normalized.  The  inteqral  over  R  is  the  integral  over  the  pulse 

volume,  and  the  V  inteqral  is  over  the  total  unambiguous  velocity  range 

of  the  radar. 

Generally  the  returned  fxawer  varies  with  position  within  the  pulse 
volume;  however,  to  incorporate  the  effect  of  imperfect  response  and  make 
the  above  relation  tractable,  a  simplifying  assumption  is  required.  It 
is  assumed  that  during  the  time  the  radar  samp'les  a  given  region  to  form 
the  sample  data  set  (used  to  form  the  Doppler  power  spectrum) ,  a  charac¬ 
teristic  distribution  (N_^e  of  particles  is  swept  past  each  point  with¬ 
in  the  pulse  volume.  In  essence,  this  states  that  the  reflectivity  is  con¬ 
stant  throughout  the  pulse  volume;  i.e., 

q(V  ,R)  -  g(V^)  (19) 

and  allows  Equation  18  to  be  re-written  as 


/2  >  ->->->■  I  /  a-  ■>->•-»  I 

V  (R) I ( r , R) dR  -I J  V  (R) I ( r , R) dR J 


where 
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This  is  essentially  the  classical  definition  of  Doppler  spectrum 

variance.  However,  the  significant  difference  here  is  that  V  (R)  is 

pD 

the  power-weighted  imperfect  response  of  the  particles,  and  thus  incorpo¬ 
rates  the  concept  that  during  the  sampling  period,  at  a  given  point  in 
space,  the  various  size  particles  composing  the  particle  distribution  ex¬ 
hibit  different  velocities  from  each  other  and  the  single-valued  air 
velocity  at  that  point.  With  an  ensemble  (denoted  by  <^  ^> )  of  observa¬ 
tions,  and  noting  that  ensemble  averaging  mathematically  is  a  summation 
and  can  be  brought  within  the  integral.  Equation  20  becomes 


<VAR>  =  <  V2  >  -  <V  >  =  <  V2  \  -  VAR- 
'  nD  x  nD  '  N  t>D  '  V 


where  now  the  overbar  indicates  averaging  over  the  pulse  volume  and  noting 
that  =  0.  This  states  that  the  average  Doppler  spectrum  variance 
equals  the  average  one-dimensional  power-weighted  particle  turbulent  var¬ 
iance  minus  the  average  variance  of  the  Doppler  mean  velocity.  This  is  a 
modification  of  the  relation  of  Rogers  and  Tripp  (1964).  It  is  important 
to  note  here  that  this  radar  mean  particle  turbulent  variance  is  not  the 
air  turbulent  variance,  but  the  power -weigh ted  response  of  the  particle 
distribution  to  the  air  turbulent  motions.  Similarly,  the  Doppler  mean 
velocity  term  also  includes  the  imperfect  response  function. 

To  incorporate  the  response  functions  it  is  useful  to  transform 

8 

Equation  22  into  wavevector  space.  Srivastava  and  Atlas  (1974)  have  shown 


Rogers,  R.R.,  and  B.R.  Tripp,  1964:  Some  radar  measurements  of  turbulence 
in  snow,  J.  Appl.  Meteor.,  1,  603-610. 

Srivastava ,  K.C. ,  and  D.  Atlas,  1974:  Effect  of  finite  radar  pulse  volume 
on  turbulence  measurements,  J.  Appl .  M-‘ cor. ,  li,  472-4HO. 


that  the  average  variance  of  the  Doppler  .spectrum  (liquation  17,  for  the 
case  of  perfect  tracers)  is  given  by 


<VAR>  =  / 4»n(K)  (1-4) 


(K)  d-4)  (K)  )dK 


(24) 


where  K  s  the  three-dimensional  turbulent  wavevector,  <?,  >  K > 


6 

=-■  ( 2 :  )  1 

4>n(K)  i 

city  com; 
the  part 
21,  Lite  1 


V  K )  F  (K)  where  K  ^  ( K )  is.  flu-  t  <»ui  i  er  f  rarisform  oi  7  (  f,bt.  , 
the  turbulent  power  density  function  of  the  radial  (1)  air  veio- 
onent.  This  relation  (Equation  24),  however,  does  not  include 
cle  response  functions.  Referring  back  to  Equations  19,  20,  and 
ew  modified  form  of  this  relation  becomes 


<VAR>  -/♦„ 
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i  the  power-weighted  imperfect  response  function  is  explicitly 
lence  R  > (K  ) ,  R  (Kz>  are  the  point  turbulent  particle  variance  and 
lean  point  turbulent  velocity  weighting  factors,  respectively, 
he  desired  result. 

c  are  a  few  [joints  to  note  in  this  relation.  First,  the  results 
ident  only  on  the  slope  (X)  of  the  particle  distribution  function 
>n  the  actual  particle  concentration  number  Nq .  Secondly,  if  the 
ume  is  decreased  in  size  to  a  jx>int,  (in  which  case  the  beam  fil- 
ion  4>j(r,R)  behaves  as  a  delta  function,  for  which  4>j(K)  becomes 
■r  all  <,)  there  exists  a  minimum  variance  which  must  always  be 
Tts  value  is  given  by 


2|K 


)  -  R  (K  ) )  dK  . 
z  1  z 


(28) 


C>(, 


Numerical  calculations  show  that  this  term  is  small  and  reflects  that 


the  two  weighting  functions  are  nearly  identical  in  net  effect.  The  exis¬ 
tence  of  this  minimum  variance  is  expected,  since  each  point  in  space  is 
observed  for  a  finite  amount  of  time,  during  which  a  distribution  of  par¬ 
ticles  having  differing  radial  velocities  (because  of  varying  response 
characteristics)  sweeps  by  and  contributes  to  the  Doppler  power  spectrum. 

The  variation  of  these  response  functions  with  wavenumber  component 
is  shown  in  Figures  E-3a  and  E-3b  for  various  particle  distribution 
slopes.  These  figures  show  that  the  power-weighted  precipitation  resjxanse 
is  good  for  turbulent  motions  having  an  apparent  scale  1  (measured  along 
the  z-direction)  larger  than  about  250  m,  with  virtually  no  response  to 
scales  of  motion  smaller  than  about  10  m.  As  expected,  better  response  is 
obtained  in  light  rain  (A  =  60  cm  S  than  in  heavy  rain  (A  =  20  cm  S 
environments.  It  should  also  be  remembered  that  the  response,  although 
excellent,  is  not  perfect  for  scales  greater  than  1  km.  Finally,  note  that 
the  differences  between  Figures  E-l  and  15—2  result  directly  from  the  diftor- 
ent  response  functions.  Equations  8  and  26. 

Equation  25  may  be  demonstrated  graphically .  For  mathematical  ease, 
the  radar  viewing  direction  is  chosen  as  vertical,  resulting  in 


<VAR>  -/<*ZZ<VVV  '♦z.f'VVV’  “z  ,29> 

where 

<}>  (K  )  =  1 1  <$>  (K)dK  dK  (30) 

zz  z  JJ  zz  x  y 
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Here  d>  (K  )  is  the  true  air  point  one-dimensional  longi tudina  power 
zz  z 

spectrum  function  along  the  z-direction.  Similarly,  ‘J)  ^(K^)  epresents 

the  spectrum  that  would  result  from  powe  spectrum  analysis  of  Doppler 
radial  velocities  obtained  from  a  set  of  successive  range  bins  along  the 
radial  (close  enough  so  that  the  pulse  volume  dimensions  may  b‘'  considered 
constant)  at  some  range  r  from  the  rad< 


67 


SCALE  X(Km) 


SCALE  X(Km) 


Figures  E-4a  and  K-4b  show  the  two  conijx>nont  Lonns  of  this  one- 

dimensional  representation  for  the  two  eases  of  outer  scale  A  =  ■»,  1  km. 

o 

The  radar  is  assumed  to  have  a  full  half-jxawer  beamwidth  of  one  deqree  and 

the  pulse  volume  depth  is  200  m.  The  outer  curve  ([>oint)  is  the  true  jxaint 

turbulent  air  spectrum  <1)  (K  )  (Equation  10).  The  lower  solid  curves 

zz  z 

(5,  20,  40,  60,  120)  are  the  corresjonding  beam- filtered  turbulent  spectra 

_  ( K  )  (Equation  31)  for  the  five  ranges  mentioned.  The  lines  » A  =  10, 
zzt 

60  cm  )  represent  the  modifications  to  the  appropriate  curve  (R^fK^) 

for  point  curve;  R, (K  )  for  range  curves)  by  introduction  of  imperfect 

1  2  -1 
tracers  in  rain  environments  corresponding  to  h  -  10  and  60  cm  .  The 

precipitation  response  weighted  curves  for  ranges  greater  than  20  km  are 

nearly  identical  to  the  respective  perfect  tracer  (air)  curves  and  are 

not  drawn. 

The  Doppler  spectrum  variance  for  the  case  of  perfect  tracers  is 
given  by  the  area  between  the  perfect  tracer  point  and  appropriate  range 
spectral  curves.  Similarly,  for  imperfect  tracers,  the  area  between  the 
imperfect  tracer  f>oint  and  appropriate  range  curves,  for  each  given  par¬ 
ticle  distribution  slope,  is  the  variance  in  real  precipitation  environ¬ 
ments.  Caution  must  be  exercised  in  interpreting  these  spectra,  since  the 
curves  represent  the  col  lapse  of  the  three-dimensional  turbulent  spectrum 
onto  the  K  axis.  However,  the  figures  clearly  show:  a)  the  loss  of 

energy  at  small  spectral  scales  (large  K  )  in  precipitation  environments, 

z 

b)  the  region  of  abrupt  change  in  slope,  the  knee  of  the  curve,  is  a 
rough  measure  of  the  outer  turbulent  scale  length,  and  c)  imperfect  particle 
response  has  only  a  small  effect  on  power  spectra  determined  from  radar 
data  whc"'o  the  radar  range  is  greater  than  about  20  km.  This  last  obser¬ 
vation  results  from  the  strong  pulse  volume  filtering  action  at  these 
ranges  reducing  the  amplitude  of  the  power  spectrum  so  much  that  further 
reduction  from  imperfect  particle  res|»nsG  is  not  significant. 

Finally,  calculations  indicate  that  in  light  rain  situations,  the  to¬ 
tal  radar  measurable  variance  is  partitioned  equally  between  x  =  y  «  z 
comjjonents .  In  heavy  rain  environments,  the  horizontal  variance  components 

are  noticeably  smaller  than  the  vertical  component  for  outer  scale  lengths 

7  2 

less  t  han  about  .26  km  (  <(  u  y  /  y  =  .95,  .91,  .  89  for  =  .5,  .25, 

. 167  km  respectively).  This  is  a  real  effect  and  results  from  the  response 
function  being  ont iroly  dependent  upon  K  .  However,  a  more  exact  function¬ 
al  'orm,  allowing  for  some  slight  direct  dependence  on  K^,K^  would  presum- 
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-  =  20  Km 
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Figure  E-4a.  Turbulent  power  spectrum  components  of  radar  Dopplei  spec¬ 
trum  variance.  Curves  labeled  pt .  represent  perfect  tra¬ 
cer  (pt.)  and  imperfect  tracer  (A  =  10  -  60  cm-*)  precipi¬ 
tation  environments.  Curves  labeled  5,20,40,60,120  km  are 
radar  beam  filtered  spot  ra  for  perfect  tracers  for  ranges 
indicated.  Full  half-|xjwer  beamwidth  is  1  degree  and  pulse 
volume  depth  is  200  m.  Outer  scale  A  1()20 


(cm5*) 


Figure  E-4b.  Turbulent  power  spectrum  components  of  radar  Doppler  spec¬ 
trum  variance.  Curves  labeled  pt.  represent  perfect  tra¬ 
cers  (pt.)  and  imperfect  tracer  (A  =  10  -  60  cm-*)  precipi¬ 
tation  environments.  Curves  labeled  5,20,40,60,120  km  are 
radar  beam  filtered  spectra  for  perfect  tracers  for  ranges 
indicated.  Full  half -power  beamwidth  is  1  degree  and  pulse 
volume  depth  is  200  m.  Outer  scale  X  =  1  km. 


ably  retain  equality  amonq  the  three  component  variances  to  shorter  scale 
values . 


The  above  results  show  that  solution  of  Equation  25  yields  the  Doppler 
spectrum  variance  that  would  be  measured  dionq  any  viewing  direction  in  any 
precipitation  environment  and  exhibiting  turbulent  outer  scale  values  of 
.5  km  or  greater. 

Figures  E-5a,  b  show  the  variance  (normalized  to  C't^  * )  versus  the 
range,  and  the  ratio  of  pulse  volume  width/depth.  The  solid,  long  dash, 
and  short  dash  curves  represent  the  air  (perfect  tracer)  and  imperfect  tra¬ 
cer  (A  =  60,  10  cm  "S  environments  respectively.  The  curve  marked  A 

o 

=  00  is  essentially  Figure  E-4a  and  is  well  represented  by  the  analytical 

9 

approximation  of  Frisch  and  Clifford  (1974).  It  is  noted  that  the  spectrum 
variance  is  essentially  independent  of  outer  scale  length  for  ranges  less 
than  about  20  km.  This  simply  reflects  that  at  short  ranges  the  various 
turbulent  fields  are  nearly  indistinguishable  to  the  relatively  small 
(compared  to  outer  scale  size)  pulse  volume  and  the  variance  contribution 
can  be  reasonably  estimated  from  the  approximation  of  Frisch  and  Clif¬ 
ford  (1974).  Significant  deviation  from  this  relation  occurs  when  the 
maximum  pulse  volume  dimension  approaches  the  outer  scale  value. 

The  variation  due  to  precipitation  environment  is  significantly 
smaller,  with  increasing  dependence  on  rain  environment  occurring  with 
decreasing  outer  scale  size  and  decreasing  range.  These  features  can 
more  easily  be  seen  by  inverting  the  figure  to  form  the  estimate  for  eddy 
dissipation  rate,  the  fundamental  turbulent  parameter. 

Figures  E-6a,  b  show  e  normalized  to  (— —■)  as  a  function  of  range, 

pulse  volume  width/depth  ratio  ,  outer  scale  A  ,  and  particle  distribution 
A.  Observe  that  the  precipitation  effects  become  less  important  as  the 
outer  scale  length  increases.  This  reflects  the  more  favorable  resfwnse 
of  the  particles  to  increasing  turbulent  scale  lengths,  which  also  contain 
increasingly  greater  portions  of  the  total  air  turbulent  energy.  Addi¬ 
tionally,  precipitation  effects  become  significant  at  short  ranges.  For 
example,  at  a  range  of  10  km,  the  assumption  of  perfect  tracers  can  lead 
to  a  32  percent  underestimate  of  e  in  heavy  rain  environments.  This  effect 

9 

Frisch, A. S.,  and  S.F.  Clifford,  1974:  A  study  of  convection  capped  by  a 
stable  layer  using  Doppler  radar  and  acoustic  echo  sounders,  d._Atiiio:;. 

Sci.,  31,  1622-1628. 
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I'iqure  1-5,1.  Doppler  spectrum  variance  (VAR)  normal  ized  to  eddy  dissi 
pat  ion  rate  (*.)  and  Kolmoqorov  constant  (C)  versus  radar 
ranqe.  R.nlar  full  half-i>ower  benmwidth  is  1  deqree  and 
inil.se  volume  depth  is  200  m.  Pulse  volume  width/lenqth 
rat  io  at  tori.  Outer  scale  A  ~  \ ,  4 ,  •*'  Km. 
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Figure  E-5b.  Doppler  spectrum  variance  (VAR)  normalized  to  eddy  dissipa 
tion  rate  (F)  and  Kolmogorov  constant  (C)  versus  radar 
range.  Radar  full  ha  1 '  -x>wer  beamwidth  is  1  degree  ami 
pulse  volume  depth  is  is.O  m.  Pulse  volume  width/length 
ratio  at  top.  Outer  scale  A  =  2,  N  km. 
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Figure  E-Oa.  Eddy  dissipation  rate  (i  )  normalized  to  Doppler  spectrum 
variance  (VAR)  and  Ko 1 moqorov  constant  (C)  versus  radar 
range.  Radar  beamwidth  is  1  decree  and  pulse  volume  depth 
is  200  m.  Ratio  of  pulse  volume  width/length  at  top. 

Outer  scale  \  =  1  <»  Wm 

C) 

76 


i 

I 


i . 
ii 


x  1000  On-') 


RANGE  (Km) 


Figure  E-6b.  Eddy  dissipation  rate  (t)  normalized  to  Doppler  spectrum 

variance  (VAR)  and  Kolmogorov  constant  (C)  versus  radar 

range.  Radar  beamwidth  ;  l  degree  and  pulse  volume 

depth  is  200  m.  Ratio  oi  pulse  volume  width/length  at 

top.  Outer  scale  1  -  .5,  4  km. 
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results  from  the  small  size  pulse  volume  being  sensitive  to  only  small 
turbulent  scales  to  which  precipitation  responds  poorly.  Thus ,  in  situa¬ 
tions  where  the  energy  source  for  turbulence  is  obtained  from  small-scale 
str  ict  ares  (e.  g.  ,  breaking  Kel vin-Helmholtz  waves)  for  which  the  outer  seal' 
may  be  small,  or  at  ranges  less  than  about  20  km,  it  is  imperative  that  th- 
pre  -ipitation  environment  be  taken  into  account. 

Of  prime  importance ,  howevei ,  is  the  thunderstorm  environment, 
era  t  observations  have  generally  found  turbulence  to  be  patchy.  Tins  is 
to  oe  expected  since  within  storms  there  exist  local  regions  of  high  shear 
which  can  provide  the  energy  for  turbulence  production.  There  is  yet  no 
way  of  measuring  by  incoherent  radar  a  turbulent  outer  scale  value  for 
these  regions.  The  usefulness  of  estimating  it  from  power  spectrum  analy¬ 
sis  of  Doppler  radial  velocity  values  is  also  questionable.  Doviak  et  al 
(1978)^  indicate  that  the  variance  of  the  Doppler  spectrum  does  not  in¬ 
crease  dramatically  with  range  for  storm  ranges  of  50  to  120  km.  This 
suggests,  after  comparison  with  Figures  E-5a,  b,  that  turbulence  may  fre¬ 
quently  be  the  major  contribution  to  Doppler  spectrum  variance  and  that, 
on  the  average,  the  turbulent  outer  scale  may  be  1  to  3  km.  This  result 
would  be  in  agreement  with  general  aircraft  observations. 

In  this  ranqe  the  variability  of  e  with  both  A  and  precipitation 

o 

environment  is  not  too  great.  It  may  be  possible  to  form  a  solution  repre¬ 
senting  an  average  A  (or  simple  set  of  solutions  dependent  on  various 

o 

estimated  A  )  with  some  simple  correction  determined  from  reflectivity  fac¬ 
tor  to  account  for  precipitation  response,  which  may  be  representative  of 
thunderstorm  environments.  Such  a  relation  would  always  be  superior  to 
the  commonly  used  relation  of  Frisch  and  Clifford,  which  as  observed  in 
Fiqurc  E-6,  will  always  underestimate  the  eddy  dissipation  rate. 

The  material  presented  so  far  deals  with  the  contribution  to  Doppler 
spectrum  variance  from  precipitation  in  a  given  turbulent  environment. 

This  variance  has  been  shown  to  be  directly  dependent  upon  the  two  basic 
turbulent  parameters.  The  eddy  dissipation  rate  is  the  quantity  of  inter¬ 
est  since  it  describes  the  turbulent  energy  available  per  unit  wavevector 
interval.  However,  Figures  E-5  and  E-6  demonstrate  that  some  knowledge  of 
the  turbulent  outer  scale  is  required  before  we  can  successfully  estimate 
this  quantity.  It  was  stated  that  turbulence  in  thunderstorms  may  be 

*°Doviak,  R.J.,  D.  Sirmuns,  D.  Ernie,  and  C.R.  Walker,  1978:  Considerations 
for  Pulse  Doppler  radar  observations  of  severe  thunderstorms,  J ■  Appl . 
Meteor.,  17,  189-205. 
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localized,  that  thei  e  are  a  tew  patches  wherein  the  turbulence  is  r*. uglily 
homoqeneous.  The  following  m, aerial  si  lows  result;.  of  ti  st.  analyses  a 
aircraft  thunderstorm  qust  velocity  data  which  appear  to  support  this 

concept  and  yield  estimates  of  both  fc:  and  A  . 

o 

The  following  data  were  acquired  on  May  26,  1976  by  an  Air  Force  F-4 

instrumented  aircraft.  Two  penetrations,  separated  by  approximately  30 

minutes,  were  made  at  a  heiqht  of  about  16  kft  through  a  portion  of  a  storm 

complex  located  40  km  south  of  Norman,  Oklahoma.  Aircraft  measurements 

were  taken  every  .1  second.  The  first  run  (Run  1)  was  from  210°.  The 

second  run  (Run  2)  passed  through  roughly  the  same  storm  region,  but  from 
o 

80  and  30  minutes  later.  Aircraft  speed  was  at  198  m/sec  during  both 
runs.  The  discussion  will  center  around  the  analysis  of  the  vertical  gust 
velocity  data. 

Figures  E-7a  -  e  and  E-8a  -  e  are  plots  of  the  vertical  qust  velocity, 
derived  gust  velocity,  environmental  temperature,  aircraft  normal  accelera¬ 
tion,  and  barometric  altitude  during  the  two  runs.  The  prominent  feature 
in  Run  1  is  the  apparent  strong  ur>ward-moving  current  between  52510  and 
52540.  It  is  about  6  km  wide  and  is  characterized  by  a  warmer  tempera¬ 
ture  than  the  environment  during  entry  and  increased  normal  acceleration 
of  the  aircraft.  Also  note  the  large  spikes  in  the  vertical  gust  velocity 
of  a  somewhat  periodic  nature,  suggesting  the  possible  presence  of  waves. 
Similarly,  Run  2  exhibits  two  aprparerit  upward-moving  ait  currents  centered 
at  53775  and  53930  sec  respectively.  The  period  53555  -  53645  exhibits 
strong  periodic  fluctuations.  Increased  aircraft  normal  acceleration  and 
spiikes  in  vertical  gust  velocity  are  also  noticed  here. 

Casual  observation  of  these  two  i>lots  suggests  that  each  gust  velocity 
time  sequence  may  be  segmented  into  a  number  of  discrete  local  turbulence 
zones.  In  Run  1,  we  focus  on  three  zones  labeled  t,  5,  and  7  in  Figure 
E-7a  corresponding  to  the  visually  estimated  turbulence  classes  moderate, 
heavy  within  upidraft  region,  and  moderate,  respectively .  Similarly,  in 
Run  2,  the  time  series  is  segmented  into  zones  1,  5,  6,  9,  7,  8,  respect¬ 
ively. 

Each  gust  velocity  time  series  was  subjected  to  pxrwer  spectral  analy¬ 
sis.  The  analyses  were  conducted  on  segments  of  the  original  data  sets. 
Each  segment  consisted  of  one  or  more  of  the  individual  zones  and  is 
labeled  1-7  (Run  1)  and  1  -  9  (Run  2)  respectively.  Each  segment  was 
linearly  detrended  before  being  subjec  d  to  spectrum  analysis. 
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Aircraft  time  aeries  data  for  Run  2:  a)  vertical  gust  velocity,  b)  derived  q  .st 
velocity.  Analyzed  segments  are  labeled  1-9. 


detrending  forces  the  power  spectral  density  to  go  to  zero  at  zero  fre¬ 
quency,  the  spectral  outer  scales  are  generally  well  beyond  the  influence 
of  this  scale  removal  action,  and  should  represent  a  true  change  of 
spectral  slope.  Figures  E-9a  -  g  and  K-lOa  -  h  show  the  |jower  spectral 
density  (PSD)  plotted  as  a  function  of  radian  frequency.  An  approximate 
space  scale  is  also  included.  Observations  from  Run  1  will  be  presented 
first. 

It  is  observed  that  the  spectra  may  be  classified  according  to  where 

the  sharp  changes  in  slope  (knee)  of  the  curves  occur.  Spectra  3,  4,  and 

6  exhibit  apparent  slope  changes  at  scales  less  than  1  km,  where 's  spectra 

1,  2,  5,  and  7  appear  to  have  the  knee  located  at  scales  greater  than  2  km. 

It  is  believed  that  segment  7  (spectrum  7)  would  also  have  shown  a  small 

outer  scale  length  if  the  time  record  52550  -  52555  had  not  been  included. 

Reference  to  Figure  E-7a  shows  that  spectra  1,  2,  and  5  are  strongly 

biased  by  inclusion  of  the  upcurrent  region  as  a  turbulent  fluctuation 

rather  than  as  a  true  feature  of  the  larger  storm  wind  structure.  The 

updraft  region  thus  acts  as  a  large,  very  energetic  eddy  which  strongly 

2 

biases  the  strength  (V  )  of  the  turbulent  field  and  introduces  extra  energy 
into  the  long  wavelength  region  of  the  PSA  curves.  In  the  sense  that  the 
upcurrent  does  not  qualify  as  a  truly  turbulent  feature,  it  is  inappropriate 
to  include  it  blithely  in  PSA.  Spectra  4  and  6  show  the  power  spectra 
derived  from  the  fluctuations  imbedded  in  the  [x>sitive  and  negative  slope 
portions  of  this  upcurrent  region.  It  is  seen  that  these  spectra  are  simi¬ 
lar  to  spectrum  3.  Referring  back  to  Figure  E-4b,  it  is  observed  that  the 
knee  of  the  spectral  curve  well  represented  the  turbulent  outer  scale  length 
under  the  ideal  conditions  of  homogeneous  isotropic  turbulence.  It  is 
felt  that  here  the  knee  may  be  a  reasonable  measure  of  this  parameter. 

Now  consider  Run  2.  These  spectra  are  not  as  easily  categorized.  It 
is  clear  that  spectrum  4,  which  includes  all  data  as  turbulence,  exhibits 
the  largest  apparent  outer  scale  value  and  generates  large  spectral  densi¬ 
ties  at  long  wavelengths .  Again,  this  results  from  inclusion  of  the  up¬ 
draft  features  as  turbulent  gusts.  The  lowei  jxawer  but  prominent  peak  near 
80  km  in  spectrum  1  results  from  a  .5  sec  ripple  riding  on  the  gust  velocity 
trace.  These  fluctuations  are  too  fine  to  be  seen  in  these  highly  reduced 
traces.  Tn  spectrum  0  a  long  (about  40  .sec)  modulation  is  likely  the  cause 
of  the  spectral  peak  near  4.8  km.  It  is  instructive  to  compare  spectra  1, 
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FREQUENCY  HI  (RRO/S) 


FREQUENCY  W!  IRRD/S) 


Figure  R-lO(a-d).  Turbulent  power  density  (m  /s  )  versus  frequency 

(r.nJ/sec) .  Equivalent  space  scale  at  top.  Labels 
(1)  -  (4)  indicate  appropriate  segments  from  Run  2 


POWER  SPECTRUM 


5,  and  <>.  Note  tliat  spectra  1,  5 ,  and  <>  have  a  wc.dl-dc  I  i  nod  knee. 
However,  when  these  three  regions  are  combined  into  one  (3) ,  the  resulting 
spectrum  does  not  exhibit  this  clear  slope  chanqe  at  short  scales.  In  this 
instance,  segment  5  appears  to  act  similarly  to  an  upeurrent  region  for  the 
composite  segment.  The  apparent  turbulent  outer  scales  for  Run  2  for  the 
individual  segments  are  iri  the  range  1  -  2.4  km.  These  data  indicate  that 
inclusion  of  storm  structure,  non-turbulent  fluctuations,  and  composites  of 
local  turbulent  patches  can  strongly  influence  the  resulting  apparent  tur¬ 
bulent  power  spectra  derived  from  the  data. 

Table  E-l  provides  the  mean  square  velocity,  estimated  outer  scale 
length,  and  estimated  eddy  dissipation  rate  for  these  various  spectra. 

Note  that  the  outer  scale  lengths  are  considerably  larger  for  Run  2  than 
for  Run  1.  This  difference  is  primarily  a  result  of  sampling  non-identical 
storm  regions,  but  may  also  indicate  a  temporal  evolution.  The  almost 

synergistic  effect  of  combining  local  turbulence  zones  into  single  larger 

2 

units  is  dramatically  demonstrated  by  the  air  motion  variance  V  .  Sur¬ 
prisingly,  the  eddy  dissipation  rates  for  the  two  spectral  sets  are  very 

similar.  Within  non-updraft  reqions,  both  runs  show  £  in  the  range  20  - 
2  3 

50  cm  /sec  .  The  largest  £  are  found  within  the  updraft  regions  and  range 
2  3 

from  417  -  1680  cm  /sec  .  These  values  indicate  that  non-updraft  regions 
would  be  characterized  as  moderate  in  strength,  while  the  updraft  regions 
are  heavy  to  extreme  (MacCready,  1964) l1 

finally,  it  is  instructive  to  consider  the  turbulent  velocity  and 
Oopplcr  spectrum  variance  estimates  which  would  have  been  made  by  a  Doppler 
radar  scanning  these  same  regions.  One  method  of  estimating  this  is  to 
apply  to  the  data  a  filter  which  represents  the  equivalent  2-way  beam  fil¬ 
ter  function  of  the  radar.  Fiquros  K-ll  and  E-12  show  the  equivalent  radar 
mean  turbulent  velocity  and  spectrum  variance  that  would  have  resulted  in 
a  uniform  reflectivity  field.  Each  gust  velocity  segment  was  detrended  be¬ 
fore  application  of  the  filter  and  the  radar  pulse  volume  was  assumed  to 
have  a  one  degree  beam  width  and  200  in  depth. 

Generally,  the  derived  spectrum  variances  are  quite  small.  Run  1 

2  2 

shows  they  are  of  the  order  of  .05  -  1  m  '/sec  with  a  mean  of  roughly 
2  2 

.1  m  /sec  in  regions  3  and  7,  while  in  the  updraft  region,  the  range  is 

* 'MacCready,  P. ,  1964:  Htandardizat ion  of  gustiness  values  from  aircraft, 

.1  ■  Appl .  Meteor.,  3,  4  39-449. 
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TABLE  E-l .  ESTIMATED  TURBULENT  PARAMETERS 


ure  E-ll.  Equivalent  radar  Doppler  spectrum  mean  and  standard  deviation  derived 
from  aircraft  gust  velocity  data.  Radar  beam  filter  has  a  1  degree 
full  half-power  beamwidth  and  200  m  pulse  volume  length.  Data  for 


gure  E-12.  Equivalent  radar  Doppler  spectrum  mean  and  standard  deviation  derived  from  air¬ 
craft  gust  velocity  data.  Radar  beam  filter  has  a  1  degree  full  half-power  beam- 
width  and  200  m  pulse  volume  depth.  Data  for  Run  2. 


m 


2  2  2  2 
significantly  larqer  at  .  1  -  8  m  /sec:  and  has  a  moan  value  of  .25  m  /sec  . 

Again,  the  spikes  are  a  result  of  the  periodic  fluctuation  in  the  original 

data  and  may  not  really  be  representative  of  turbulent  contribution. 

The  results  for  the  late  period  are  similar.  Zones  1  and  6  have  var- 

2  2 

iances  of  the  order  .01  m  /sec  with  the  updraft  regions  exhibiting  .25  - 

2  2 
8  m  /sec  . 

These  equivalent  radar  variance  values  are  generally  smaller  than 

2  2 

those  typically  observed  by  radar,  which  are  generally  1  -  36  m  /sec"  with 

2  2 

a  mean  value  range  of  4  -  16  m  /sec  .  The  discrepancy  here  may  result 
from  the  vertical  qust  velocity  not  being  representative  of  the  horizontal 
turbulent  components,  and  exclusion  of  normal  radar  spectrum  broadening 
agents  such  as  system  noise  and  wind  shear.  The  preliminary  analysis  in¬ 
dicates  that  it  is  really  necessary  to  compare  coordinated  aircraft  and 
radar  data  to  determine  whether  equivalent  radar  estimates  can  be  re¬ 
liably  obtiiined  in  this  manner. 

The  results  discussed  here  will  allow  for  increased  accuracy  in  Doppler 
radar  estimation  of  turbulence  intensity,  and  will  have  practical  appli¬ 
cation.  As  earlier  stated,  the  two  main  contributions  to  Doppler  spectrum 
variance  are  the  transverse  shear  of  the  radial  wind  and  turbulence.  The 
Doppler  radar  can  estimate  the  shear  component  by  observing  the  variation 
in  azimuth  and  elevation  of  the  mean  Doppler  velocity  between  successive 
radials.  Thus,  the  turbulent  contribution  can  be  estimated.  With  some 
knowledge  of  the  turbulent  outer  scale,  and  storm  reflectivity  factor,  a 
reasonable  estimate  of  the  eddy  dissipation  rate  may  be  obtained.  This 
information  can  then  be  combined  with  the  known  response  characteristics 
of  aircraft  for  identification  of  regions  of  hazard  to  aircraft.  Use  of 
these  results  witli  incoherent  radar  is  significantly  more  difficult,  how¬ 
ever.  The  incoherent  radar  can  not  estimate  the  Doppler  mean  velocity,  and 
therefore  will  not  allow  lor  the  determination  of  the  transverse  shear. 

Thus,  the  turbulent  variance  contribution  may  not  be  well  known. 

12 

Atlas  and  Srivastava  (1971)  proposed  a  method  for  estimating  an 
atmospheric  structure  function  with  incoherent  radar  by  forming  combinations 
of  signal  returns  from  successive  pulse  volumes  along  a  given  radial.  They 
showed  that  the  mean  squared  difference  of  Doppler  velocity  between  the 

12 

Atlas,  D.,  and  R.C.  Srivastava,  1971:  A  method  for  radar  turbulence  de¬ 
tection,  I EKE  Trans,  on  Aerospace  and  Electronics,  A ES - 7 ,  No.  1,  179-187. 
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adjacent  pulse  volumes  can  be  related  to  the  DoppJer  spectral  variances  by 


2  A2  (P1  +  P2) 

(v2  "  V  =  T  ~Tjr—  <VAK12  <pv  +  'V  '  VAK'i  '  i  "  VAH,  ,'2)  (  n) 

where  A  is  the  radar  wavelenqth,  VAR  and  VAR?  are  the  Doppler  spectral 

-2 

variances  (sec  )  trom  pulse  volume  numbers  1  and  2  having  total  spectrum 
powers  P^  and  P^,  and  VAR^  is  the  variance  of  a  coni|>osite  Doppler  spectrum 
obtained  from  adding  the  siynals  from  the  individual  pulse  volumes.  Un¬ 
fortunately,  this  relation  is  a  highly  smoothed  estimate  of  the  true  atmo¬ 
spheric  structure  function  (Sychra,  1974)'*"*  requiring  very  careful  inter¬ 
pretation;  it  is  also  related  to  the  environmental  wind  shear. 

A  measured  Doppler  radial  estimate  may  be  considered  as  a  sum  of  a 

large-scale  background  wind  (V  ),  turbulent  com|«nent  (V  ),  and  an  addi- 

w  t 

tional  random  error  and  noise  bias  (e) ,  as 


v  =  V  +  V  +  e. 
w  t 


(34) 


With  this  convention  the  left-hand  side  of  Equation  33  becomes 


<V2  -  V*  ‘  <V2»  -  Vl»)2  *  <»»  -  »,/  *  ‘“2  - 


(35) 


where  the  second  term  is  actually  the  true  Doppler  estimate  of  the  pulse 
volume  filtered  turbulent  structure  function.  Tf  two  azimuthal  pairs  of 
pulse  volumes  and  (V3,V2))  arG  combined  in  this  manner,  then  in 

a  region  where  the  mean  turbulent  air  velocity  over  the  region  occupied  by 
the  pulse  volumes  varies  little, 


(v3  -  V  -  (v2  - 
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Using  Equation  36,  a  linear  variation  in  J  urge-sea ) e  wind  can  bo 
fitted  to  the  data,  thereby  providinq  an  estimate  of  the  large-scale  shout. 
If  a  quadratic  variation  in  the  wind  field  is  desired,  two  such  sets  of 
Equation  36  may  be  used. 

Obviously  this  technique  is  prone  to  error  where  the  qradient  of 

13  -  ■ 

Sychra,  J.  ,  1974:  Fluctuation  .spectra  and  velocity  structure  function. 

Tech .  Rpt . _ 3_3,  Lab.  for  Atmos.  I’rob.  ,  t  of  Chicago. 
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moan  turbulent  velocity  is  lame.  However,  at  the  r.inqcs  of  interest,  per  - 
hups  60  200  km,  most  ot  the  turbulent  energy  is  mapped  into  tile  Do:.!  lor 

spectrum  variance.  For  this  case,  the  radar  estimate  of  the  turbulent 
structure  function  is  small,  and  its  variation  is  not  significant,  allow¬ 
ing  Equation  3f>  to  be  useful.  The  validity  of  this  relation  is  independent 
of  the  actual  magnitude  of  the  air  turbulent  intensity.  Therefore,  the 
above  re  ations  may  bo  applicable  within  local  patches  of  turbulence  and 
may  return  reasonable  estimates  of  transverse  radial  wind  shear,  thus  per¬ 
mitting  an  estimate  of  eddy  dissipation  rate  to  be  obtained.  Usefulness 
of  this  technique  can  only  be  established  after  comparison  of  radial  shear 
estimates  from  Doppler  and  incoherent  techniques  using  real  data. 

(C) .  Conclusions 

The  results  of  theoretical  investigation  demonstrate  that  useful  and 
understandable  relationships  exist  between  Doppler  radar  spectrum  variance 
and  turbulence.  For  ranges  less  than  20  km,  Doppler  radar  spectrum  var¬ 
iance  is  independent  of  turbulent,  outer  scale  and  dependent  upon  eddy  dissi¬ 
pation  rate  and  range.  At  ranges  qreater  than  20  km,  Doppler  spectrum 
variance  is  strongly  dependent  upon  all  three  (quantities.  However,  once 
the  lamest  pulse  volume  dimension  becomes  comparable  to  the  turbulent 
outer  scale  length,  the  spectrum  variance  becomes  range  independent.  Im¬ 
perfect  particle  response  is  found  to  reduce  significantly  the  measurable 
turbulent  strength  for  outer  scale  lengths  less  than  about  500  m.  At 
larger  lengths,  precipitation  resqxjnse  has  only  a  minor  effect.  It  is 
suggested  that  order  of  magnitude  estimates  of  turbulent  eddy  dissipation 
rate  may  be  obtained  by  using  a  generic  outer  scale  length  of  2  km. 

Accurate  estimation  of  turbulent  strength,  however,  awaits  some  method  of 
remotely  estimating  turbulent  outer  scale  values  for  the  local  turbulent 
patches  within  thunderstorms. 
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